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THE PHONON BOTTLENECK IN OPTICALLY PUMPED RUBY 
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ABSTRACT 


The resonant interaction between radiation and two- 
level states, with application to the interaction between 
29-cm + phonons and excited cr3t ions in optically pumped 
ruby, is studied via the Heisenberg equations of motion for 
the system. As a first example the spontaneous decay of a 
Single two-level state is studied. Next a set of rate 
equations, rigorously valid for low densities of two-level 
states, 1s derived from the Heisenberg equations. The rate 
equations show that in the case of far-from-equilibrium 
initial conditions the system relaxes rapidly to an inter- 
mediate quasi-equilibrium with the density of electronic 
States in the upper level increasing with the density N* of 


neo : 
two-level states, 1.e. cr? 10ns in the case of optically 


*% 


pumped ruby, like N tomelargerNe wee lterc shown that 
Spatial diffusion is ineffective as an escape mechanism for 
the phonons for large N*. Thus the subsequent decay of the 
system is controlled by spectral diffusion of the phonons 
away from the resonance frequency due to repeated absorp- 
tion and emission by the two-level states, with accompany- 
ing rearrangement of the phonon and electronic energies in 
the coupled mode. The phonon spectrum inside the excited 
region of the ruby crystal (with diameter b) is shown to 
suffer self-reversal and to broaden with time to a final 


iS 
width proportional to (N*b)’*. The effects of the inhomo- 


geneous broadening of the resonance, due to possible 
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random static strains in the ruby crystal, on the evolu- 
tion of the electron-phonon system are also briefly 
examined. Finally some exact solutions of the Heisenberg 
equations are found and improvements on the rate equations 


suggested by field theory are studied. 
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CHAPTER 1 
INTRODUCTION 


1.1 General Survey 

The interaction of radiation in resonance with a sys- 
tem of two-level atoms has long been a subject of interest. 
Resonance radiation refers to that part of the spectrum 
that is emitted by an atom in an electronic transition from 
some excited state to the ground state or a metastable 
state. If most of the atoms of the system are in the ground 
state (or the metastable state) then this radiation is 
highly absorbable and will very likely be absorbed by 
another atom after traversing a short distance, thus excit- 
ing that atom. This situation is referred to as the 
"trapping' of the resonance radiation. By a ‘random walk' 
process of emission and subsequent absorption by different 
atoms the radiation eventually reaches the boundary of the 
system and is lost into the bath. 

If a phonon rather than a photon is emitted in the 
electronic transition, the resonance radiation trapping is 
referred to as the phonon bottleneck. In this case the 
decay of the two-level atoms leads to a transfer of energy 
to a narrow band of phonon modes near the resonance fre- 
quency Vu, which in turn deliver this energy to the thermal 


Datheoteciersunroungimoscoystaie selLiethestatewoL decay of 
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the atoms is slower than the rate of delivery of energy to 
the bath by the phonons, then the phonon modes near reson- 
ance act aS a reservoir i.e. they are almost always empty 
and can absorb all the energy emitted by the decaying atoms. 
We may say that the phonon modes are at the bath tempera- 
ture of 0 K. Thus the excited two-level atoms can decay 
unhindered with their spontaneous decay rate 1/T,. On the 
Other hand, if the rate of delivery of energy to the bath 
by the phonons, due to multiple absorption and emission, is 
Slower than the rate of decay of the atoms, then the phonon 
modes near resonance will become effectively full. Before 
another atom can decay, it must wait for some of the phonon 
modes near resonance to empty their energy to the bath. 
Because of this phonon bottleneck the effective decay time 
ae of the system is prolonged. 

At first glance it might appear that since phonons 
are bosons an unlimited number can be put in any mode so 
that the two-level atoms should always be able to decay 
freely. However statistical considerations show that the 
phonons and atoms reach a quasi-equilibrium with the frac- 
tion of phonon modes near resonance that are occupied never 
exceeding the fraction of two-level atoms that are excited. 
We may say that the phonon modes near resonance, as mea- 
sured by the occupation numbers, heat up to the temperature 
of the two-level atoms and no higher. 

The earliest experiments on the trapping of resonance 


radiation were done on the optical resonance lines of 
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mercury and sodium vapor under low pressure. For example 


hy 
Webb in 1924 and Zemansky> Liev Ze soe LeAmLiICe trapping Of 


Sar Dey ON photon emitted in the ip Soa ground state 


transition of the mercury atom. They observed that the 


trapped radiation decayed in times of the order T. #10 sec 


(which is Gig times longer than the lifetime T, ~10 sec 


of the 6°, state), and that T. increased with the density 
Nano mercury atoms in the ground state for small N* as one 
would expect. However at the largest densities N* they 
LOUnds that T actually decreased very slowly with N*. This 
rather surprising result was attributed either to pressure 
broadening or to the more frequent collisions at higher 
pressures cauSing non-radiative transitions from the radiat- 
ing 6°P, 


lower. In the former case, collisions with the radiating 


state to the metastable 6°P stacesur2USe volts 


atom cut off the radiating wavetrain, thus increasing the 
natural linewidth (hence the name pressure broadening) and 
effectively reducing the lifetime T, of the excited state 
and presumably also the decay time T Of the system. In 
the latter case, a competing process simply reduces the 
number of atoms in the radiating 6°P, State. Aside from 
these factors one must take into account the hyperfine 
structure of the 2537 A line, the Doppler broadening of the 
line due to thermal motion of the gas molecules, and a host 
of emission and absorption processes induced by atomic 


collisions. Thus we see that optical experiments on gases 


are very 'unclean' i.e. there are many complicating factors 


that make interpretation of the experimental results diffi- 
cult. We shall see shortly that many of these problems can 
be eliminated by doing similar experiments on magnetic 
impurity doped crystals at low temperatures. 

The first theoretical treatments of resonance radia- 
EPLOnmLLAaApOLNOmWweLemclVenm by ak. ls Compton” ine! O22 eandsb. 
bia in 1926 in connection with the transmission of light 
in the atmospheres of stars. They assumed that the radia- 
tion quanta propagated from atom to atom by a type of 
random walk. Thus the density of excited atoms N, could 


2 


be described by the diffusion equation 


Z 


ON, (r,t) ' 
—. & V N, (x,t) ’ (1.1) 


ot 
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where A is the mean-free-path of the quanta, and 
Clea 


is the average velocity of the quanta taking into account 
the time of the collision i.e. the excited state lifetime 


Ty as well as the propagation time between collisions sain 


In the case of photons one can obviously neglect a in 


comparison to T The treatment of Compton and Milne how- 


L’ 
ever did not’ take into account the fact that the mean-free- 


path of the quanta is not constant but highly frequency 


-= 
' 


7: abet Z t Agta sie , ifew 2) 


devendent. The mean-free-path of radiation on resonance is 
very short whereas that of radiation far from resonance is 
very long and limited only by the dimension b of the con- 
Cato Dele emmcaninguotu\inehdee) 3)e1S, unclear. 

bre Pod 7 Pole tetog derived an integral transport 
equation describing the time decay of the density of 
excited atoms N, (r,t) taking into account the frequency 


dependence of the mean-free-path, namely 


aNA(c,t) Na (x,t) 


~ Ont 
0 


Q 


The first term on the right hand side of (1.4) describes 
the spontaneous decay of excited atoms at position r in 
time Th: The second term describes the excitation of atoms 
at position r due to radiation coming from other atoms cy 
These atoms 7 decay in time Ty 


spectrum f(v). We assume that f(v) is strongly peaked 


emitting radiation with 


about the resonance frequency v and normalized to unity. 


roca 
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The radiation is attenuated by a factor e 
to absorptions in travelling from r' to r (here a(v) =I1/A(v) 
is the frequency dependent absorption coefficient), and by 
a factor 1/4n|r-r' |? due to spherical spreading of the 
photons outward from ‘ee The radiation reaching position 
ris then absorbed with probability Cia Fee nal LV yewealntLe— 
grate over all atoms rc’ and all radiation frequencies v. 
Implicit in Eq. (1.4) is the assumption that the radiation 


travels from r' to r instantaneously 1.e. that Wai Ue 
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Thus this equation is appropriate for photon but not phonon 
radiation trapping. 

Honeeeaine solved Eq. (1.4) by assuming a solution 
exponentially decaying in time and using a variational 
method to find the associated time constant. Vek lenko/ in 
1959 provided an alternative solution of Eq. (1.4) which 
made use of the so-called Ambartsamyan transformation. [In 
both cases the analysis is very complicated and the results 
cannot easily be generalized to arbitrary emission spectra 
eA 

We wish now to present a novel and very simple solu- 
tion of Eq. (1.4). We can transform the integral equation 
(1.4) into a differential equation® by assuming that No is 


not too rapidly varying in space so that we may make a 


Taylor expansion of N, (r',t) about No (x,t). Substituting 


N,(x',t)*N y Wuker.) (ri-r,) x 


2 deh 
1j 


dil 
Wes — 
Bisbee (2a r,) V,Nj (r,t) +5 


a N, (r,t) (Sie 25.) 


under the integral in (1.4), we immediately get the result 


@oebe 
ON. (r,t) ws 
Tia Ceca Ga N, (r,t) ’ (1.6) 
where 
i 
a=] shy sE( ONO) 5 (ia) 
0 Bie (v)T) 


which is simply the frequency dependent diffusion coefficient 


1/307 (v)T averaged over the emission spectrum f(v). Eq. 


i 
(1.6) is just the diffusion equation (1.1) with a now well- 
defined effective diffusion constant. Neglecting factors 


2 


Of order one we may replace V' by Lyne where b is the 


dquametecrao Let nescontalner andetind the result that N, 


decays with time constant 


& EY) 
T. = “foe : ((aPyais9) 


(a (v)b) 


As soon as the emission spectrum f(v) is specified (1.8) 
cam Dekinteqratedror @o(Vjielsiproportional to £(v) din®ther- 
mal Soar yan accordinomctopkherchoLtes Lawl} mEorgexample 
Toen()iyistagdeltagtunction dvev 5) or a rectangular line- 
shape, then 


2 
Hee T, (a (vo) b) ’ (Le) 


which is the usual result for material particles. For 
Doppler and dispersion lineshapes we find that 


+ 
Prism smericetoceamletreatenctes y Lor which ay) > 1/b° 


On the other hand for those frequencies v far from resonance 
ELOuUsWwiLChePLOpOLULOnal ty would amply that (Vv )e<.1/D,awe 
must set a(v) =1/b since for those frequencies the mean- 
free-path A(v) =1/a(v) is limited by the size b of the 
COntaIne Gael ill SmeloenOLm I neCcOntLadloulOnmsOtm it TChOtt Bou law 
because these frequencies do not come into thermal equili- 


brium due to the absence of collisions. 


dads T, a (v5) b Yin (a(v.)b) Clee SO") 
and 
T. ~ T) vatv.) b (eee 1s) 


respectively, in agreement with the results of Holstein’ 
and Veklenko’. 

Since we will derive an expression similar to (1.8) in 
Chapter 5sewe deter a detatled discussion of it until then. 
For now we simply wish to remark that Holstein and Veklenko 
CrdsniOumebyeLOmGi1scOVer = Lie physical basis of the different 
Seed melicomeleiopliyeancdethateras 61.4) could not give 
them any information on the radiation spectrum inside the 
container or emerging through the boundary of the container. 
We will try to answer these questions in Chapter 5. 

A completely different semi-classical formulation of 
the trapping problem was given by Jacobsen and Stevens” 
(1963) who were interested in the effects of a crystal 
lattice of paramagnetic spin 1/2 atoms immersed in a mag- 
neti CmrLeld= On thee propagatrom of ~phonons “inthe crystal. 
Their system is described by the Hamiltonian (we set f=1) 


2 
P 2 
A= }{lA 871+ [sn+S (ne Udy ines (Umea ee dcalian(Ual2) 


2 
n 
where the sum over n is over’ all lattice sites in the one- 
dimensional crystal; the first, second, and third terms in 


the square brackets are respectively the spin, phonon and 


: 1 
Spin-phonon interaction energies; Pw Une and Sn are 


Mespect lve yeticumomenctum, displacement, and eee component 
of spin of the tie atom; and A , K and g are respectively 
the Zeeman splitting of the spin energy levels in the 
external field, the nearest neighbor restoring force, and 
the coupling between the strain and the spin at position n. 
The Heisenberg equations of motion for Pw Us and 32 
(Similar to the ones we derive in Chapter 2) can be com- 
bined to give a pair of wave equations for the displacements 


and spins, namely 


ee Zz xX eS Ok 
mU Se rR 2U,) = tose Sane . (oi 3a) 
"xX 2s 9) 1 x z, 

a ae ae pens = “a e e 
neRrs sit OA Been gA (U4 OSes Cia} 


a 


Here T, =T, (9) is the phenomenologically introduced life- 
time of the spin against phonon emission. With no inter- 


action (i.e. g=90 and hence =o} (l.l3a) 1s the wave 


EL 
u 
equation for the propagation of free phonons and (1.13b) is 
the equation for an undamped harmonic oscillator of fre- 
quency A . With the interaction switched on (g and T, 
finvee) = thesrigcht wand side Ot Bdqs..(l.l5a,bjeact -asescOources 
anGmsinkse tom tNespnononssandespins aE lecmicerachlOlmdcs 
, ; Z 
non-linear due) mathematica la yetomthetiactor ela Le Ae 
and due physically to the existence of processes such as 
stimulated emission, etc. ‘To solve Eqs. (1.13) we replace 


om by the c-number average <S7> Cen OGean Orla a DOOLa= 


tion (i.e. few excited atoms) or acs Se / 2 LOL an 


inverted population (few unexcited atoms), and assume that 


U, and Ss" bothwvarvelike et Wt-kna) 


where na is the posi- 
tion along the chain of the ae atom. This yields the 


dispersion relation 
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is the shifted resonance frequency, and c, =avK/m and 
c(w) =w/k are respectively the speed of sound (phase 
velocity) for zero and non-zero spin-phonon coupling. In 
Figure 1.1 we have plotted the real and imaginary parts of 
co /c (a) =d6(w) -ia(w), corresponding respectively to the 
dispersion and attenuation of phonons of frequency o. 

In Fig. l.la the solid curve representing a normal popu- 
lation (<S*><0) has negative slope in the region around 
A' indicating anomalous dispersion where a group velocity 
cannot be defined. Outside this region the increasing 
positive slope aswapproaches A' indicates a decreasing 
group velocity of the coupled spin-phonon wave. This is 


due to the delay by a time T, of phonons every time they are 


il 
captured by normal atoms, the resonant phonons being cap- 
tured most often. On the other hand, for inverted popula- 


tions the dashed curve of Fig. 1.la indicates an increase 


in the group velocity near the resonance frequency beyond 


trae 
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the phase velocity ome This however is in contradiction with 
a result by Sommer fe1a?® who shows that it is not possible 
to propagate signals at velocities greater than Co: Figs 
1l.lb shows that for normal populations (solid line) the 
attenuation a(w) of the phonons is a maximum near reson- 
ance, whereas for inverted populations (dashed line) there 
is actually an amplification of the phonons near resonance. 
This is completely analogous to amplification by stimulated 
emission in the laser and has been demonstrated experimen- 
Bed dewey Tucker! + ri Of 

Thus far the formulation of Jacobsen and Stevens 
describes the attenuation and dispersion of incoming plane 
waves by a system of two-level atoms. However it is un- 
clear how to proceed and describe the evolution of the 
two level atom-phonon system from non-equilibrium condi- 
tions because of the anomalous dispersion near resonance 
for which we cannot define a physically meaningful group 
velocity. Moreover their formulation supposes that the 
paramagnetic ions form a regular lattice whereas we are 
interested in the situation where the ions constitute 
impurities randomly distributed throughout the crystal. 

There are other rate equation models of the phonon 


Bocoene cheat wa 


, but virtually without exception they 
assume that the emission f(v) is rectangular. This is an 
unrealistic’ assumption and from Eqs “1(1.9)— (Cl e11)” we see 


that the decay time T. of the system depends crucially on 


the lineshape. 


1.2 The Phonon Bottleneck in Ruby 

In our discussion of the experiments on trapping of 
optical resonance radiation in gases we mentioned that the 
experimental results were difficult to interpret because 
Of a variety of effects such as collisions of gas atoms, 
etc. Many of these effects can be eliminated by studying 
instead the trapping of resonant phonons by paramagnetic 
TORNMIMPUGI tL 1eSsinecrystals at low Pempers clire sagas aabecsuce 
the ions are then not free to move and collide. We will 
be mainly interested in ruby, which is an A1,0, crystal 
doped with cr ion impurities, but our discussion is 


applicable to many other impurity-doped crystals such as 


Sapphire which is Al.0, doped with a ions. 


oeo 
The relevant part of the electronic energy level dia- 
gram for the cr?* ion in ruby is sketched in Dies eee The 


free ion ground state is a with L=3 and S=3/2. The seven- 


fold orbital degeneracy is split in a cubic crystalline 


Ere Ldn tO-a “A, orbital singlet and the ti. and in trip- 


lets which are shown cross-hatched and are the main absorp- 
tion bands of ruby. The famous fluorescent or state arises 
from the 2G level and lies below the aT and 0, levels. 

The presence of a small trigonal crystalline field in com- 
bination with spin-orbit coupling further splits the i, 
level by 0.4 cm + and the 2a state into two so-called 
Kramers doublets 2A and E separated by approximately 29 cme. 


The 2A and E states are metastable against optical transi- 


tions to the ground state with lifetimes of approximately 


ue 


Axial field+t+ 
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Electronic energy levels of the Cr “ 


Taken from Reference 19. 
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4 msec for both. The 2A state also decays to the E state 
by the emission of a 29 cmt phonon in a lifetime Ty ofa 
few nanoseconds. In the phonon bottleneck experiments of 
interest to us the 2A and the metastable E levels serve as 
the upper and lower levels of the two level atom system and 
then29 cm + phonons serve as the resonant quanta. In 
Figure 1.2 we have also shown that each doublet is split in 
an external magnetic field. 

Let us now describe some of the experiments which 
have studied the bottlenecking of the 29 cm phonons 
resonant between the 2A and E levels of the cr?* DOn@1 ii 
ruby. The most complete experiments were done by K.F. Renk 
and G. Bani 4 ames ca. 1977. In their experiments a cw-Argon 
laser at 514 nm is used to optically pump a small cylindri- 
cal region with diameter b=.7+1mm (full width of the 
laser beam at half intensity) and 1 cm long (length of 
crystal) inside a l cm? ruby crystal (doping .03 wt3% cr3*) 
immersed in liquid Helium at 2K. This causes the excita- 
ELonEof emer ions to the “Ty and “2, levels, and by radia- 
tionless decay, the population of the 2A and E levels in 
approximately equal numbers. However in the steady state 
established by the continuous action of the laser nearly 
all of the excited er3t ions are in the E state because the 
lifetime Ty of the 2A state is much shorter than the life- 
time of the metastable E state and the thermal population 


Otethes2Ametatcmtclowsati slow crystal temperature. The 


emission of Ry and Ry Lactation ma (Cimelid.. see etcan, bemused 


Ls) 


aS a weak probe for the measurement of the population of 
the E and 2A states because the lifetimes of these states 
against optical transitions to the ground state are nearly 


equal and so long compared to T The cylindrical excited 


L° 
volume containing the crt ions in the E state is addi- 
tionally pumped with the radiation of a pulsed dye laser at 
580 nm with a pulse length of 2-3 ns, a repetition rate of 
30 Hz, and a power of 50 kW. The schematic arrangement of 
the experiment is shown in Fig. 1.3a. This pulse leads to 
an initial non-equilibrium population of the 2A level. It 
leads as well to a time-dependent population of the E level 


which however is small compared to the permanent population 


N* of the E level, as evidenced by the fact that the Ri 


Signal remains more or less unchanged, and that the R, Sig- 


l Sronal. = The Ri and R, 


Signals are observed with nanosecond time resolution using 


nal is small compared to the R 


several time-to-pulse-height converters in parallel. From 


the ratio of the time dependent R, signal to the time inde- 


Z 


pendent R, signal one can obtain the relative electronic 


V 


population NSC) Neer: Ce) / Re (where N, and N*tare the) den— 


sities per cm? Oe cr?* ions in the 2A and E states respec-— 
tively). It is important to note that the laser pulse is 


weak enough that the condition N. << N® is always satisfied. 


AAV 


In their experiments Pauli and Renk measured the 


intensity at the Ri and Ry lines following the laser pulse 
for various concentrations n* (i.e. for various intensities 


of the cw Argon laser) keeping constant the total initial 
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(a) Experimental arrangement for high densities Osa ene. 
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Fig. 1.3. Schematic experimental arrangement used by Pauli 
and Reng ee fOmRENCeS cudveO te chem phonon bottle— 
neck in ruby. All signals entering and leaving the 


ruby crystal are transferred via light pipes. 
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energy injected into the electron-phonon system! (i.e. the 
dye laser pulse energy). 

In Fig. 1.4 we show typical decay curves that they 
LOundsLorethe R, Signal. The results shown are the counts 
accumulated after approximately 10° pulses. In the inset 
we show the initial fast time behavior. We see that after 


the dye laser pulse reaches the sample, the 2A levels 


peak 
SD 


the R, Signal. Phonon emission sets in within the lifetime 


become populated instantaneously to a peak value in 


Ty Seetioe2Nectatesleadung to a tastidecay of Ry 


time of approximately 3 ns (the duration of the laser 


(ey) aan fi 


pulse). Ry (t) does not vanish but rather falls to a value 
Boe because many of the phonons do not escape but are re- 


absorbed in resonance by other er?* VOnsSminethesfestate . As 
a result a quasi-stationary state of the coupled electron- 
phonon system is established, which we will henceforth 
refer to as the quasi-equilibrium. The electrons and 
phonons are not in complete equilibrium, but evidently 
always close to equilibrium as the coupled system slowly 
decays away. Some of the phonons decay anharmonically or 
escape from the excitation region and are quickly replen- 
ished in 2A+E transitions as indicated by the slow decay 


inethe:keestona Weinwos Lime T of the order of microseconds. 


2 


rn this thesis the term ‘electron system' refers to the 


system of ert ions in the E state (in the cylindrical exci- 
tation region) which can be excited to the 2A level by ab- 


SOLDELOnGOL san29 em 2 phonon. 
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The R5 panes aSsa, Lunction, oc seimestound .by 

Pauli etwal at different concentrations N* of 
excited ae ions in the E state. In the inset, 
an Ro Signal BtatiNey 25 Saskia bmn) as plotted 


linearly on a nanosecond time scale. 
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From the R, Signal one can extract the relative popu- 
lation No tN" of the 2A level after the fast decay and the 
decay time T of the slow decay. These quantities are 
plomcedmine rig eal Sauand bas funce@ons of the density N” 


of E levels. One finds that 


nen constant for N~ < 10 cm , 
Ne Nee ae ty D5epe—3 to) 
nN* 2 for Nae LO cm ; 
and that 
constant for N* < nde cm > ’ Cl. 17) 
a C= +h 1) 5) . 
net? fOoreNas>) 100. cm ; 


Evidently the phonon bottleneck first appears at a density 
Né= 10> ai. At lower densities the phonons are not 
trapped and escape more or less freely with the ballistic 
flight time while at higher densities the phonons are 
trapped for increasingly longer periods of time. 

For completeness it should be pointed out that for 
low densities N” < 1073 oye Renk and Pauli used a heat 
pulse to inject the initial energy into the electron-phonon 
system (seer big. U.3b) rather than a dye laser pulse.) ) ihe 
heat pulse was produced by passing a 100 ns pulse of elec- 
ECUCECUGLeltstnbOugimcameLe a) mane film) of Constantan evapor— 
ated on the surface of the ruby crystal. It was estimated 
that a Planckian distribution of phonons is produced (with 


Pemperatucers 2Ush saErOUCi Lyena leo ieWhd Cheenter the ruby. 


Crystal and half of which are lost to the®helium bath. 
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Fig. loo (a). The) quasi-equi librium population N5I/N ps eka: (es) 


the slow decay time 1h ee eats: 2A states as functions 


of the concentration N* (ene Cage ions in the E state 


20,21) The circles were obtained 


frome opticalwexcitation ssa Seal the squares 


found by Pauli and Renk 


from heat pulse enaarexetharnmieer > omits excitation region 


diameters b=.7 mm and b=1.9 mm respectively. 


ee. 


The advantage of the heat pulse method is its high repeti- 
tion rate (104 pulses per second vs. 30 pulses per second 
for the dye laser). It is, however, unsuitable at higher 
densities N* > or em > because the resonant phonons are 
not able to penetrate well into the interior of the excita- 
tion region due to resonance scattering. 

In another similar experiment Pauli and Renk util- 
ized a single cw-Argon laser to set up a steady-state popu- 
lation of excited cr ions and phonons and found that No7 
No t/N™ increased linearly with the excitation region dia- 
meter b. This indicates that the phonons escape spatially 
and not by anharmonic decay as had been previously sugges- 
ted. In the anharmonic decay hypothesis the resonant 
phonons are so effectively trappnved that they decay in time 
ooh into two or more non-resonant phonons which can then 
escape ballistically. Orbach and Vredevoe“* and Rlenensss 
have done theoretical calculations which indicate that 
‘eget ns and 380 ns respectively for longitudinal and 
transverse phonons. However Kaplyanskii et ihe aie in time- 
of-flight experiments found that the phonons propagated at 
least 1 cm with little attenuation indicating that 
on >> 1 usec. Thus we conclude that we may neglect 
anharmonic decay of phonons in our analysis. 

Lengfellner et aie 22 studied the lineshape of the 
2A+E transition by splitting the 2A and E levels in an 


external magnetic field as shown in Fig. 1.2. In their 


experiments optical pumping at low crystal temperature 
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results in the population of the Ely state. An infrared 
laser at 29's7 cm + (slightly greater energy than the 
difference A between the 2A and E states in zero magnetic 
field) is used to illuminate the crystal. As the external 
magnetic field is increased from zero, the 2A and E levels 
are both split (see Fig. 1.2) by increasing amounts and at 
two values of the magnetic field Ry radiation is detected 
indicating athat Es 2A and then E , > 2A transitions 


a! 1 i * 
2 +74 2 


nr 


have taken place. The R, Signal as a function of the mag- 


Neciceurretd Nasa. Lorentz—like lineform for both transi— 


Pronsmnatcatingml1 fetimesbroadenmng) of the 2A level) (half 


di 


width Av = 012 sat) and hence a lifetime T = (27Av ) = 


L 
-44 ns. However Rives and Meltzer?? have pointed out that 
there may be a Significant contribution to the linewidth 
from inhomogeneous prosdening: implying a longer lifetime. 
Geschwind et Ai? also performed an experiment which 
measured the lifetime T) indirectly by observing the transi- 


At low temperatures the direct decay is very 


. =x ae _— 
eon Ei E 
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slow and the transition takes place in two steps via the 


T Random local strains in the ruby crystal cause locally 
different crystalline electric fields and hence locally 
different splitting between the 2A and E levels. Macros- 
copically one observes a Gaussian distribution of Lorent- 
zians i.e. a Voigt profile. This smearing is known as 


inhomogeneous broadening. 
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and 2A.,-+E (The 


so-called Orbach process: E +h +h 3 


+, 


Raman process, in which the intermediate state is virtual 
and may lie at any distance above the ground state below 
the Debye frequency is a second-order process and not impor- 
tant at these low temperatures). A theoretical analysis 


shows that the relaxation time for the E level is given 


il 
+35 


by 


flip A/kT 
Torbach ~ e 2 


L 
enh 


(igels ) 


where A=29 cmt, T is the crystal temperature and eee 


is the time for the spin-flip transition 2A, 75 2. con— 
co 7 
parison with experiment yields nee 


ner 


=15 ns. Theory pre- 


dicts that the non-flip transition is 50 times faster or 


Tana y= 23 ns so that the combined decay rate of the 
DMelevel is ight, YE a ee? aan sees et) Stee te 


should be mentioned that Kurnit et al. in phonon echo 


oer menesa: found that T, = 2-8 ns and Rives and wetenee= 


in time-resolved experiments similar to those of Pauli and 


Reni ca ae previously described found that T, =1.1 ns. Thus 


we conclude that there is some uncertainty in the value of 


Pegi tcl incul .Otethc 2A state. 


To conclude our survey of phonon bottleneck experi- 
ments in ruby we briefly mention the results of two other 
groups. Rives and elie as already mentioned, 
performed time-resolved experiments similar to those of 
Pauli et alz and) found qualitatively similar results. 
However rather than measure the concentration N” of ey 


ions in the E state they used N* as an adjustable parameter. 


- ‘ tae ~ ene 


Tilswltetssotrriculteto make a quantitative comparison of 
their experiments with theory. Dijkhuis et edb 
studied the steady state properties of the bottleneck 
using continuous pumping to populate the E and 2A levels. 
They found (see Fig. 1.6a) that the intensity R, varied 
quadratically with Ry indicating strong bottlenecking of 
the optically generated 29 em)” phonons. (We shall show 
later that a linear dependence would indicate no bottle- 
necking.) Furthermore they found that when an external 


Magnetic field is turned on the R, intensity and thus the 


2 
bottleneck is reduced fourfold as shown in Fig. 1.6b. 
Evidently this is due to the 2A+E transition being grad- 
ually split into 4 distinct transitions by the Zeeman 
effect (see Fig. 1.2). The number of phonon modes to be 
bottlenecked thereby increases by a factor of 4 and the 
bottleneck decreases by a factor of 4. From the Gaussian 
shape of the curves in Fig. 1.6b Dijkhuis et al. conclude 
that the 2A>E transition has a Gaussian lineshape, possibly 
due to inhomogeneous broadening. Note that the lineshape 
is the broader the larger Ri isestNatwethiseindicates that 
the phonon spectrum is broadened at large N47 se oThe insetnos 
Fig. 1.6b gives the external field 6H at which R,/R, (H=0) = 
0.5 (this is a measure of the width of the phonon spectrum) 
vs. Ri: 
This concludes the general introduction and survey 


of experiments on radiation trapping and the phonon bottle- 


neck in ruby. 
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10 0 (b) Magnetic field (kG) 


Ca)) Intensity of R line (counts/s) 


Pie. Ona) he intensity of the R, line vs the intensity of 


line found by Dijkhuis et a imecontinuous 


ene Ri 
pumping experiments. The slope “of 2 indicates ‘strong 
bottlenecking of optically generated 29 cmt phonons. 


(Dye Then intensity "of the R. line, normalized to zero 
eS comialmmagneticmrer1ecld, vVSylthe magnetic field ator 
several excited state populations N" as measured by 
the R, intensity. The inset gives the field SH at 


a 


which Ro/R, (H=0) =0.5 vs Ri: Both figures taken from 


Reference 14. 
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1.3 Aim and Scope of the Present Investigation 

The objective of this investigation is to present a 
theoretical analysis and interpretation of the phonon 
bottleneck experiments in ruby; notably those done by G. 


Pauli et a caus and Dijkhuis et Nile ees 


which we briefly 
aeScribeduinwulea:. 

We begin in Chapter 2 by studying the electron- 
phonon interaction and deriving a field theoretical 
Hamiltonian describing it. After a survey of the solutions 
to the interaction problem which have already been studied 
in various approximations, we close the chapter by deriving 
a set of linear equations of motion for a set of quantum 
operators describing the system. 

In Chapter 3 we study the simplest case of a single 
atom interacting with a phonon field. We study in detail 
the spontaneous decay of the single two-level atom as an 
example of the decay of an unstable quantum system. 

In Chapter 4 we make the transition from quantum 
field equations to classical rate equations in parallel 
with W. Pauli's famous derivation of the master equation. 

Chapter 5 constitutes the bulk of the thesis in which 
we solve the rate equations and give a detailed comparison 
with the results of the bottleneck experiments. 

In Chapter 6 we try to improve the rate equations by 
going back to the quantum equations and extending them and 


by presenting some exact solutions of them. 


CHAPTER 2 


bole DD THEORETIC DESCRIPTION OF 


THE PHONON BOTTLENECK (T) 


We wish to present a rigorous theoretical analysis of 
the phonon bottleneck in ruby, at each step of our argument 
clearly stating the assumptions and approximations necess- 
ary to arrive at our conclusions. Basic to the understand- 
ing of the phonon bottleneck is an understanding of the 
electron (or two-level atom)-phonon interaction, which we 


now discuss. 


2.1 The Electron-Phonon Interaction 

Let us begin by deriving a Hamiltonian for our system. 
As far as the ere ions are concerned we are interested 
only in electronic transitions between the 2A and E states. 
Transitions between these levels and the ground state or 
the absorption bands are used only to monitor or to feed 
Chessy stemeana dO nNotrakrect the dynamics of the bottie— 
neck, so they are ignored. Thus we may approximately 
describe the Qn ion at position rasa two-level system 
and use the vectors ae and ae to represent the ion in 
the 2AgandiEistates respectively. Let us introduce the 


Pauli spin-1/2 operators 


om 2 i 0 0 ik ede 
sta= fesio“en | . es ml mL) 
2 0 0 ae 0} = 0 -1 
ig 1g ig 
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which we define to act (only)on the gee ion, with commuta- 


tion relations 


ats Peete nee Z 
ae et 2S. oe v 
ae Gay 
_— Zs pS = 
Sp Speci = S. See 
Then S represents the excitation transition of the eee 
ion EAL, namely 
0 oi 
+ 
ad Lm i 
and similarly Se represents the decay 2A+E 
ele 0 
S = “ Ge 4) 
5 i 1 
ig 
The levels 2A and E are separated in energy by the amount 
A = 29 om > Shy, where i rehalte GHz. If we define the zero 


of energy to lie halfway between the 2A and E levels then 


As” is the energy operator for the ieee On; asinee 


~ 


l ib 
AsZ | = 1 4| | (2.5) 
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i.e. the eigenvalues of Ase are the energies of the upper 


+ 
and lower levels. Thus the Hamiltonian for the free er? 


Fons “is 


Bb) 


gah. =} As ; (ons 
ie 


where the sum is over all excited Creuiions in the excita- 
tion region. 

Let us now turn to the lattice. We wish to analyze 
the vibrational energy of the lattice in terms of the 


excitation of phonons. The crystal symmetry of ruby is 


rhombohedral 7? with its unit cell containing two molecules 


of Al,0,. Thus if we perform the usual normal mode analy- 
Sis on the lattice vibrations we find dispersion curves 


A Similar to the ones sketched in Fig. 2.la. 


~ 


Thevexactmshape of thercurves and widthsaof- the .Brivllouin 
zone will depend on the direction k that we choose. For 


any given direction k there are 3x 10=30 branches due to 


(L) Coa) 
ie mT ’ 


~ ~ 


the 10 atoms per unit cell. Of these, three (uw 


Tr 
and 2 in Fig. 2.la denoting the longitudinal and two 


~ 


transverse branches) are acoustic branches linear ink for 


f0™) (Om) 
ae siti hy 


~ ~ 


branches. We approximate the dispersion curves using the 


small k. The rest, denoted w ares optical 
Debye model for the acoustic branches and the Einstein 
model for the optical branches as shown in Fig. 2.1b. The 
Brillouin zone is replaced by a sphere of equal volume so 


that 1 tsenadausmLs Ky where 
3 

De 

2 
T 


Sie (2.8) 
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and nis) the density of unit cells per com? in the ruby 
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crystal. The acoustic branches are replaced by straight 
lines with slopes (i.e. longitudinal and transverse sound 


VerOcmiCLe ci) = Cen C =r and C the same as the original 
L Ty T, 


curves near k=0, and the optical branches by lines with 


zero Slope (i.e. w), =constant). The average sound velocity 


c in the Debye model is defined by the relation 


Spee 
ns] = mE] + = + nd oro!) 
L Ty) T, 
Where Ce=7HC=s, mandac are averaged over all directions in 
L Th Ts, 
Space because the Debye model is isotropic. For ruby 


Cha Fag. Oe m/sec and Ky zs iG meerieae so that the Debye cutoff 


frequency is =ck, * 10x A, where Ae ¢ rn is the reson- 


D 


ance frequency. Since, as we see in Fig. 2.l1b, all the 


Optical modes have frequencies Ww, > w Cn yVeacousticeand 


k Da 


no optical modes can be excited in a 2A>E transition. 


Thus we need only consider the acoustic modes and write 


the free Hamiltonian for the phonons as 


a: 
k,|k|<k, Stalglealat = ¥ : 
| (s)" 
where k, = 2am, /L pele ei Ome with), m, an integer. by and 


Bae! are the usual phonon creation and annihilation opera- 


tors satisfying the commutation relations 


(s) (S.4) 

i Pic It an °kk! vet é (2v 1615) 
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Let us now study the Grau ion-phonon interaction. In 
the ruby crystal a small number of nae ions are substitu- 
ted by cr?* ions. Their nearest neighbors are 6 025 tons 
which form a distorted octahedron’® around the cr?t Lon, 
giving rise to a cubic crystalline electric field and a 
weaker trigonal field. These static fields, whose poten- 
tial we denote by Vor uncon UuncCtlLongwith the spin-orbit 
coupling give rise to the 29 cm + Bolitting,of the 2Avand 
E levels shown in Bie eel NmAccuL. Ol cOnchis Static 
field, vibrations of the lattice, through the charges of 
the ligands (the om ions), produce a modulating effect on 
the crystalline field and therefore on the energy levels of 
the cr3* ions. This causes electronic transitions in the 


ous ion. The total crystal field potential can be 


expressed in powers of the strain Bye 


SSP We wae apeiat. G2 7aL2)) 


ee Tea. 5 


Vorystal 


where ve as discussed above is already included in Hon! 
Ee = SOME Meg) WES is the direction-averaged strain due to 
lattice vibrations at the site of the Cras ion (U(x) is the 
crystal displacement from equilibrium at position x) and 
Vi and V5 are coupling parameters which depend on the level 
of excitation of the cr?t Lon. 

Let us consider a simple model of a monatomic one- 
dimensional lattice with sites indexed by the integer n 


with a few Col ion impurities at sites n=r. Combining 


(2.7) and (2.10) (specialized to a one-dimensional lattice) 


rita “7 on rh i nm o@ei 
22 -AOSIUEASTIAL Ter 2=rithe ] Hy: Vass we & | 


mn 


autitedrs SIE |nGr “Sr” bey gee Piene A. lates ¢Au5 end 
, 44 & A a : 


<—_> £ : 
“aa are gsodtulen feemeen, alee Sie. ~ go Se 


<= 


aoidw 


; on hee 
etlgtod Patciosgse.l EB taaee 


> > - ‘ c eo 7 fr = 4 - re 
fieli 22 . oni bierngse 220 6 2 


r 
be ~~ s + qa) 
a uy “ > a 
7% Ce eal a atte 
7 * 
L 
} A ‘ | é d 
. - i) oft Bric r 
A - . ; 
* = — £ al _ 
; os bs" s 
© 4 
. ~ era >. < ™ P « 
rt ‘ q : - re | i.e 4 Ti 
Ps 
“ { 
"ies 
« + 
. "i j> = sf r L aytee 
: r 2 = - 
: L 
> oe ‘ = . - - be 
NS 
== = 
5 >, 96. 4 Bin 
= ‘ 
~ E en ~ * - é ~ ‘.- 4 
= ea: a Ge > Ws a? w ; 74 = dee bl ‘ 


: - — ~ , ; = - - be 
ar Bi Cx) vi) et >) | S¥_8 Bibs & ema a8 
- “ ; 


bam. (y gota bNeey OS mur tilings mot +f apse ohae Ms’ {ese bl 


B dedi eek eekiayes Sa a, ks Ve 
_— ft har cera es 
anal “a0 9a er > a ) 


34 


with the first two. terms of (2.12) we get the Hamiltonian 


= A on Soh * a ; eee) 
where 

Hion =) 4s d 

is 
# = J w bib (2.14) 
ph t Kw K , ; 

and 

Vic = 1) (Cee RUSE Tee aoe 
Here ¥ iS a coupling constant, one ay is proportional 
£0. thes strain, at-the ae lattice site and V is for now 


Wie 


unspecified. In terms of the one-dimensional phonon anni- 


hilation operator 


3 See ) eK (p - imo 


- Un), (2.15) 
k (2mw,,N) * n + 


and its hermitean conjugate b) (Po and UL are the quantized 


momentum and displacement of the fe lattice site), we can 
write 
2 Sige Ke +ikr 
TUCO Tt pect Wheat oF p g,, (bye + ble lan; (216) 
where 
L 
1 = 
She = aaa, sink (57a) i ei WB 
k 


The simplest form of V that can cause transitions of the 


lex 
Chie . ws ; 
G ion and is hermitian is 


Vi —— Siuate > oka S a (2.18) 


With this form we finally arrive at the Hamiltonian for 


the electron-phonon system in one dimension 


= Zz ¥ "Tee aK Yo tel iee 
LES jp Wee, |) nites sp teams) bee (eos, < BORS Touma! 
Me k oak 
- ikr oe eed kr 
+ {Sb e + b Se iy (2-19') 


This Hamiltonian can be trivially extended to the general 
case of a three-dimensional lattice by including a summation 


over phonon branches i.e. by everywhere replacing bey ) 
Ss 


k 
ee 1’ Toi thus allowing for different 


coupling to each branch. However since experiment does 


and 3.79 s=L, T 


not give us these couplings we may as well simply choose 


them all equal. 


2.2 Survey of Solutions to the Electron-Phonon Interaction 


Problem 

THesianmyetOntLane (2. L9)eissa generalization of a 
Hamiltonian that was first studied in 1954 by Dicke? ihe 
connection with the spontaneous emission of coherent radia- 
Lion, pl. ec...» Superradiance! ,Winsaylaser. | Since: then al very 


Varges body of iaiteratuneshas Aenea 


studying the 
Hamiltonian under a variety of conditions in various appro- 
ximations and uSing various methods. 


An approximation that is almost universally made is 


the so-called rotating approximation in which 
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the second term in braces in (2.19), the anti-resonance 


Germ, sOLOpOLtlOnda.:, CO 
Sieomc os Wey Tat C200} 
iG 10, 


DoeaLoppeCa.= INemisudl=juStIriCatiron is that it does. not 
lead to first order processes that conserve energy (or 
number of excitations). To see why it is called the anti- 
resonance term and why it is negligible, we use the inter- 


action picture where 


i 0, (t) = (0, (t), Hod ’ 
and CLANS. 


i|p,(t)> = A’ |p, (t)> : 


where 0. and vy are any operator and state vector in the 


interaction picture and oe coe reonon and ap Ug 


Recor noe tou(22ijeecnemresonance terms of vary, like 


b ~i(A-o, )t i ~ 
by (t) S(t) = e b, (0) S_(0) r 
C72. 2.) 
+i (A-w, )t 
ot oe k i 
S(t) b, (t) = e S (0) b, (0) ; 
and the anti-resonance terms like 
G -i (Ato, )t » 
S(t) b, (t) =e S_-(0)b, (0) E 
C25 23)) 
+i (Atw, )t 
+ i ral k + 4 
S(t) b, (t) =e S (0) b, (0) A 


pince Grand) Os are both defined positive, for w, ~A the 


k k 
resonance terms have an approximately stationary phase, 
whereas the anti-resonance terms have a rapidly varying 
phase ee lhusewhensweerntegqrate® (2.22), fore A’ not too 
strong, we get a strongly resonant denominator for the 
former (p< L/(o-A) } and anti-resonant denominator for 
the latter (Wo 1/ (aw, +A) ) which we can neglect in com- 
parison except for large times. 


Sik 


Swain has studied the case of a single two-level 


atom interacting with a single photon (or phonon) mode of 


frequency w, with #' including the anti-resonance term. 


n 


He found that if initially the atom is excited and there 


is no phonon present, then the probability N,(t) that the 


4 


atom is excited at time t is given in terms of a continued 


fraction, namely 


Ray (2. l Bis oe tat 2 
a 2 2 
Or Key 
wa, - 29° 


ar, G2me2o 


SE 


so that the atom never decays; and that keeping only the 


first approximant to the continued fraction in (2.24) gives 
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| 1 f 3 eet | 


’ (220) 
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which results in a cosinusoidal N, (t) Vasms Own ei srl io 3 0 
OSe Clap rer) eee vceoLudy 3(2. 26) Binwmore detail in Chapter 3 
where we find thet 1t 1s the exact solution to the one atom-— 
one phonon mode problem with no anti-resonance interaction. 
Let  ussinterpret, thesegresults. .Eq.. (2.25) shows. that. the 
uncoupled atom behaves like a harmonic oscillator of fre- 
Guch Cyn COUDLINGstncwatom to another harmonic oscillator 


of frequency w, (the phonon mode), causes a shift in the 


1 
frequencies of the two oscillators (determined by the poles 
of the integrand of (2.26) ) and causes the energy of the 
coupled oscillators to be transferred back and forth cosin- 
usoidally between them at the beat frequency (= g) as for 
al eOalit Ome CcOUDleOsCclassucaleharmonic oscillators. Finally 
Eq. (2.24) shows that the inclusion of the anti-resonance 
term for small times is negligible for small coupling g but 
that for large times it causes N, (t) to lose its perfect 
periodicity because of the infinite number of poles in the 
integrand, of (2.24.4 . 

Beyond the almost universally used rotating wave 
approximation which we just discussed, the approximations 
that can be made depend on the conditions of the system. 
For example when applying the Hamiltonian (2.19) to the 
laser, usually only one radiation mode on resonance is 


considered. Furthermore if the wavelength of the radiation 


is long compared to the dimension of the container of the 
two-level atoms, when we may also drop the phase factors 


nee 


ea - F : : . 
e meas and e in (2.19) and get the Hamiltonian studied 


by Diickea ; namely 


H= bb + J AS’ +g 


- + 
(b oe + S,.b) . C2 aL) 
Us r 


Dicke considered all N atoms as a single quantum system and 
applied to it the rules for the addition of angular momenta. 


Let us define the operator 


102 
Ren ae ee CR) (2.28) 
Lex py yz 
where 
a a as (2.29) 
Ie 


is the ee component of the 'spin' of the total atom system. 
Then ne has eigenvalues r(r+l) and the operator R” has 


eigenvalues m satisfying 
i} 
im] 3 = @ 5k (2530) 


r and m are integral or half-integral according to whether 


NeVvSMeVeNROL OCG. LE N, and Ny are the number of excited 


and unexcited atoms respectively, then 
Ne=oN eee N (273 1) 


2'73 2:) 
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and=m as called ene coordination number. Note that the 


energy of the atom system is 


ao 


oe 
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Di at ga ae >) : G25) 


A straightforward algorithm using raising and lowering 
Operators can be used to generate the eigenstates [ie hus of 
Z 2 5 : ; 

R andR . Then taking matrix elements of the interaction 


energy one finds that 


<xv,m| W'|r,mtl>¢o f(r F m) (x + mtl) G2es4) 
are the only non-vanishing ones. Taking the squares of 


the matrix elements with the lower sign gives the spon- 


taneous decay probabilities 
Testo trem) (o- mt ly (2200 


where Ty is some constant. These results are conveniently 
tabulated in Fig. 2.2, where we show the states eaM> 
their degeneracies and their energies (the height in the 
diagram). Let us interpret these results. By setting 
N=1 and r=m=1/2 in (2.35), it is evident that I, is the 
BadtatlOnmerd temo n amsing veuexci ted atom. Lf ee Ne i.e. 


2 


if all N atoms are initially excited, then 


bes RN lias, C2nenoOF) 


which is just N times the rate for a single atom. The pro- 
portionality to N indicates incoherent emission by the 


acOnS eels r=5N and m=0, i.e. if half the atoms are 


excited, then 


Lai 
T= 5N(5Nt 1) 1, G2F ous) 


Al 


r=4 r=N_j r=p-2 os 
N . Bee = ae aaa 
aepre Ne 2) Ae a 
NI N-l fold degenerate Z v4 
N 
=m TD te al a Sa ae degenerate 


3(N-2)T 2(N-3)I 
fe) fo) 
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Fig. 2.2. Energy level diagram of an N-atom gas in the Dicke 
model, each atom having two energy levels separ- 


ated by energy A. Spontaneous radiation rates are 
N 
=i) 


indicated. The atom system energy is EB, A(mt 5 


Taken £Lrom Reference 30. 
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Wh1lChe GomeLarge, Neils proportional to Ne! This indicates 
coherent spontaneous emission. Dicke calls this the 
"Superradiant' state. 

Radiation trapping emerges naturally in this form- 
alism. Suppose that exactly one atom of the gas is 
Hiselailyeexclced we Ihnen the gas radiates initially at 
the normal incoherent rate but after a short time ceases, 
with the probability of the photon being emitted being only 
WNeamel is fol lowsittrom Figker? 12:0) Thecassumed <initiad 
state with m=-5 +1 is an equal superposition of the single 
state with r=N/2 and the N-l degenerate states with r= 
sol. The atten N-l states cannot decay and ‘thus the pro- 
Dabisimity got) hadlatuonet rappingweils <. (This is exactly 
analogous to the famous case of helium in which inter- 
combination between the triplet system (orthohelium) and 
the singlet system (parahelium) is not allowed by the 
selection rules so that e.g. there is no yg Pes Grane 
tion). By field theoretic methods in Chapter 6 we will 
EinGgmthatwethemarapoing result just found is exact. 

Let us now discuss some other work done on the 
electron-phonon interaction. Schart-- studied the Dicke 
model and found that the mean photon number showed oscil- 
latory behavior given by elliptic functions. Prakash and 


Ghandvane wrote the Dicke Hamiltonian in the so-called 


Boson representation by replacing 
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- iF 
Rie a sa, 
+ ‘4 
Ra A544 
2c) 
Z ly st oi 
Rv 7 (a5a5 - a,a,) 
| 
eres = IN, Nn, > 
where n, =3r-m and n, =xr+m are the number of excitations, 
and ay and a, are boson operators, which annihilate an exci- 


tation in the lower and upper levels respectively. Simpli- 
fications result if we have a highly excited assembly of 


atoms because then n, is large and the a, or upper mode 


Z 2 
Cane bet Greateds Classically. “Similariy Vt ny is large the 
ay mode can be treated classically and if n, * Ns Poe vd 


(Superradiative condition) then both modes can be treated 
classically. 

Weulibeay made computer studies of the evolution of 
the density operator of a system containing a small number 
of atoms and found oscillations in the photon field similar 


to those of sche urest 


Bek also studied various aspects 


Bonatacic ‘etlal. 
of the Dicke laser model. 

All the authors mentioned so far assumed the presence 
of only a single radiation mode and neglected the spatial 
dependence of the system. As we know from the results of 
Holstein (Eqs. (1.9)-(1.11)), the existence of a phonon 


band and the size of the system is crucial to a quantita- 


tive description of the phonon bottleneck. Also, in light 


of our discussion after Eq. (2.26), it seems likely that 
the oscillations in the models of Scharf and Mallory are 
artifacts of the single mode assumption. 


EVES! studied the Hamiltonian 


Leonardi and Persico 
Gite OLmd ACODSenn and Stevens” in connection with the 
relaxation of an inverted system of paramagnetic spins 
interacting with a phonon field. However, as we already 
mentioned in 81.2, this model assumes that there is a spin 
On each lattice site and therefore is not an appropriate 
description of the phonon bottleneck in ruby. 

Finally we should note that a number of authors have 
derived and solved macroscopic and kinetic equations des- 
cribing systems obeying the Dicke Hamiltonian (2.27) or the 
Bi vietamitecontans (2.19) For example eerste faa derived a 
kinetic equation for the photon density matrix of a monomode 
laser within the framework of the Prigogine-Résibois theory 
of non-equilibrium statistical mechanics. Grinberg and 
iceman and Bukhbinder et se derived kinetic equa- 
tions for the relaxation of a paramagnetic spin system 
using Zubarev's method of the nonequilibrium eR hater) 
operator. And D'Yakonov and Beene, using a Feynmann 
diagram perturbation method, derived kinetic equations for 
thes time evolutionsofethne density matrixgoft ansexcited acom 
when there is diffusion of resonance radiation. In case 
the density matrix is diagonal, their equation reduces to 


thatwote Hols ted hen(ho miei: 
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In this survey we have seen several qualitative fea- 
tures of the phonon bottleneck. Examples were Dicke's 
radiation trapping due to the selection rule against inter- 
combination of states of different total ‘spin' r, and 
Jacobsen and Steven's anomalous dispersion of phonons near 
resonance. The only quantitative theory however is that of 
Hows ce rosand ecomomy appropriate: coOsphoton ditfusion. 
Thus in the next section we begin to derive a quantitative 


description of the phonon bottleneck from first principles. 


2.3 Linear Equations of Motion for the Electron-Phonon 
System 
Let us begin by writing down the Hamiltonian (2.19) 


(we set h=1), namely 
eb ibe + (lige casera hte. <>) (2.39) 
+L Oy ok ) L I k e ae PK : : 


In this equation the summation over k implicitly includes a 
Summation over the 3 types of acoustic modes. Using the 
commutation relations (2.2) and (2.11) we can immediately 


write down the Heisenberg equations of motion for the oper- 


ators Dee Swrand Se namely 
hiss ws < - ike 

iby, Dd Py + d Se 7 

1 Sep aye ST ao) Gees el kt (2.40) 
c a Leth ym Karek 
2 Z vA Ry yeeetes + ikr + - -ikr 

is) = [SH] =I I, (S bye piceaDy See yee 
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It is easy to show that the number of excitations in the 


system is constant, i.e. } bb 5) ag 
in 
grating Eqs. (2.40) formally we get 


= constant. Inte- 


=f, 8 ~iker | thay, (GEAEM) 
lt) Se aay priate | at enor. S(t ae 
ie = 0 
t 
k-r 
a eats ak ee — TNC =the ; 
om) =e SEO | EST | dt'e 08 este ’ 
~ 0 
= a q , t 
Silt) = 82(0) ~4 J a4 (¢ [ at’st (tty, (e4 
= 0 i os 


~ik-x f 
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’ iy ' = ' 
she Fygils WSeAas eke (22415) 


0 


These equations can be solved by iteration. 
Let us first find the solution in first order pertur- 


bation theory. Keeping terms up to order ae we find that 


= ees) = G2) z (3°) 
SG) Sa (i) Sea L are Ge) by ge KAR, K (t) 
7 (4) o2zo- i (5) dat 
x Bethel tanta che i Scale MoM Aye es Spby Pies a 
(OB 4D) 


where all operators on the right hand side are evaluated at 


time t=0 and pe yet) (eel) et, Dea vemscOomearunc tLonseof 
time. The probability No (t) that the oe atom is excited 


~ 


at time t is given by 
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No (t) = oy SMe se (t) |i (2.43) 


~ ~ ~ 


Wierem iy micmticrstatesot the syctemmat time b=—0. After 


some algebra we find that to order gi? 


4° Lae af ph + ph 
ay AS) aie hg k F, (t) ) (Ne Chea 0) NO eee Lg 
7 (2.44) 
where 
Dit) = <ilp! (t)p) (4) |is (2.45) 
k k ik 
is the number of phonons in mode k at time t and 
A-w® 
Sin (—x—~ t) 
Saag amaee 


We wish to point out several features of the perturbation solu- 


BuO OLESiIa lst imes = lmerer( tt) = a so that if there are 


Ones 
no phonons initially present, the excited atom decays like 


Nee tyme NCO laa eng (2.47) 
k 


We note that initially the atom does not decay exponen- 


tially; rather it begins with zero slope. Indeed in the 
case of a single phonon mode the factor iega ue igcy Gelats. 


beginning of the series expansion for the cosinusoidal 


evolution of the excited atom found by Swain??. Notice 


also that meas) = is independent of k i.e. any phonon 


~ 


regardless of its energy can be emitted by the atom 
initially in accord with the energy-time uncertainty 


relation AEAt > m. However as time progresses F(t) 


~ 
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becomes more peaked around oO. =A. Leteus now assume. that 


~ 


the phonon modes form a continuum so that we can replace 


| H 
~1 


>. | DiD)-dG™, (2.48) 


where D(w) dw is the density of phonon modes per cm? in the 
frequency interval dw about w and 2 is the volume of the 


crystal. Then in the limit of very large times 


lim =F AS Oe Sy: (2.49) 
{70 k k 


expressing the fact that for large times energy is con- 
Seuved wer Luggung (2.48, 49) “into (2.44) we find that the 


probability that the atom is excited at time t is 


Ne (t) =, (0) + reine (0)nh" (0) - WN, (0) (1 +nBP(Q)) 2 (esas 
where 
oe 27 g% D(A) 9 (051) 


is the spontaneous decay rate in agreement with Fermi's 
golden rule. The three terms in braces represent the 
excitation by initially present phonons of the atom in 
the lower level,spontaneous decay, and stimulated decay 
by phonons of the atom in the upper level respectively. 

It is interesting to note that in thermal equilibrium 
where all occupation numbers must be constant, Eq. (2.50) 
yields 
(2'52) 
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Combining this with the fact that at temperature T, the 


number of atoms in the + and ¥ states are given by N, = 


4 
exp (-E,/kT) and Ny = exp(-E//kT), with 1 foe hs A, we get 
the Planck distribution for the phonons, namely 
ph 1 
“4 2 G2) 
A Qh/kT a 


Perturbation theory has the shortcoming that, although 
it gives the decay rate of an excited atom, it does not give 
directly the spectrum of phonons that is emitted. For 
exanple, nq. (2.00 )mclaims™ that only a phonon of energy 


exactly wo, =A can interact with a two-level atom. But we 


k 
KnOwmcLOMecnemunceLrtainty relacion wnat actually ali band of 
phonons of width Aw » va = 20 g* D(A) Cane interact. with 
i 


the atom. For this reason we wish to iterate Eqs. (2.41) 
cCOnalivorders. # ToOmtacilitate this we represent the equa- 
elon omoleLOG all wrOmeOw p,srOLger the phase factors, 
summations and time integrations, and focus attention on 
ClesDOWeImOLe Ne sCcOUDLIingmeg, and the order of the non= 


commuting operators. We write symbollically 
ey) x fey) Ae Gi (te) 
Smit) Su COEREEgGa (tu) bit) (2.54) 


ae 


Sate Sa Omg (Sittin CEU OM (cl) .Sm (tn). 


which we represent pictorially as 


b(t) = A C2 55a) 


and 
Z a 
Say ies E (Zoo C) 
Here the dots represent operators at time t=0 e.g. 
Die 
O = baie=07 , (2.56) 


open circles represent operators at finite times e.g. 


Or = VeGe ah ls, C2) 


vertical lines represent addition, e.g. 


= Fie Sah] : (2.58) 
Pe We aa 
and horizontal lines represent multiplication, e.g. 
aint) 2 el 2G 
‘jp Seen reg: Love ae th, (2e097) 
ae eed — teres Boel = 
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The diagrams are read from top to bottom (in ascending 
POwers Of gg). The order from left to right is important 
Since the operators in general do not commute. To iterate 
BQS.5 (404l))picctorially we simply plug the diagram for 
operator q into et, For example iterating S (t) several 


times we get the perturbation tree shown if Fig. 2.3. 


The open circles at the ends of the branches indicate that 


we have to iterate ad infinitim. To regain the first order 


perturbation result we cut off the 'tree' at the qe level 
and multiply out the branches to get 


S(t) as BA eEAa5 = -1G ialy IS See, ie (Oa 


to order g 
Re-inserting all the summations and time dependencies we 


agalimoetsexacGryenG. L024 2). 


Linearizing Approximation 


Let us now make the important assumption that there 


are very few excitations present in the system, i.e. 


}<byb, > + }<s1s/> ae esecke x ) 1 26 15) 
k ve 16 1s 


Where -<...> =<1!2..ia> 5 j%&Undeed, leesusvassume that there 
is just a single excitation initially present in the-sys- 
tem. Then it is easy to see that we may replace the per- 
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i.e. we everywhere replace aie) SS because from the 
other terms of S7(t) (see Eq. (2.55c)) there will come con- 
ee UC OM SascO S(t) [is with at least two annihilation 
operators acting on a state |i> (with only one excitation 
present), and these will vanish. Furthermore any term 

of the series for Sa) represented by (2.62) acting on 


state |i> has the form 


z . 
Se Gro ce CO}BD CO) a> 


= I ~n as 


or (22563) 


SG) sae to WS. . CO eae 


which either have the value (a2 ee or vanish depending 
on whether or not state |i> has a phonon present in mode k 
or the ype atom excited. We may summarize this by saying 

thateutethe sstatemiaemhace ustwamsinglemexcttation present, 
then we may everywhere replace Ay eile We expect this 
to be approximately valid also for a small number of exci- 


tations present. With this replacement, Eqs. (2.40) become 


S| 


. =— bok «i 
Si. Saar oF L he sy 
ie (eno) 
oa +ik.r 
ae = AS, = p Fe PRE cas : 


Whereas Eqs. (2.40) were non-linear in the number of opera- 
tors, these equations are linear and thus exclude processes 


such as stimulated emission. 
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CHAPTER 3 


SPONTANEOUS DECAY OF A SINGLE TWO-LEVEL ATOM AND 


THE GENERAL THEORY OF THE DECAY OF UNSTABLE QUANTUM SYSTEMS 


Let us first study the simple case of the interaction 
of a single atom with a phonon field. In this case Eqs. 


(2.64) reduce to 


sek = Pk a Ty 


(She)! 


+ — 
PiSGemeneCwOPpCEotOEse. and) Se refer tosthe, single atom. 
These equations can be uncoupled by the use of (one-sided) 
Houriermtranstforms.) DefinesthemFourier transforms of.S~(t) 


and by (t) to be 


§ (x) = ¢[S (t)] = | at s hey) eee ' 
0 
and (32) 
= = oct 
By (x) = fib, (t)] = | atp, (t) : 
- as 0 ae 
with inverse transforms 
Ses: Al (Ne 3 4 ee) Se 
6 (3573) 
Seal eee St 
b, (t) = F [B, (x) ] = -7E $ dxB, (x) © ; 
a; a > 


where the integration contour 6 encloses the lower half- 


plane as shown: 
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cereals 
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x-plane 
G 
Corao) 
PReCeEGULLer cranscrorms OL Eqs., (Gel), are 
te ) = ake 
ee a tae Shs (UO Usaaisc i [Oa (3.4a) 
i So = 
Shee 
3° (x) ===>, S (0) + ) —= B (x) (3.4b) 
x-A K Te : : 
Using (3.4a), dg (x) can be written in closed form as 
we it ( ta j ib; (0) 
19) (.0:) se g a : USre>y) 
ge | BS ee See 
ee Jes i . 
k k 
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Expanding the denominator in a binomial series about x-A 
gives the same series that would result from the simple 
iteration of Eqs. (3.4a,b). The atom-phonon system may be 
viewed as a single oscillator of frequency A (the atom) 
coupled to a set of N harmonic oscillators of frequencies 


Oy (the phonons). The coupling shifts these N+l frequencies 


~ 


slightly to new vaues given by the positions of the poles 
of 3 (x). There are two quite different cases depending 
on whether the phonon modes are discrete or form a 


conkinuum. 


3 lee DOCS aOR NWO nCve laa COM Lome kena temcnysita L 


In a crystal of finite extent it is well known that 


there are a finite number of discrete phonon modes. It is 


: wr 
a Oo) Seb 0 J 
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easy to show in this case that S (t) and hence the proba- 
Duepty N, (t) of the atom being excited at time t has some 
sort of periodic behavior. To be specific assume that at 
time t=0 the atom is excited and there are no phonons 

present. Dropping therefore the second term in braces, we 
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It is well known that the N+l eigenvalues of the symmetric 


matrix M, which we denote Wir Drees, aremrea linus 


t 
N+1 
inverting the Fourier transform we find that the probability 
N, (t) that the atom is excited at time t is given by 


= Pe thi eye 
MIME) = Lee hel Rey a aS } (3.8) 
k=1 


which is obviously periodic (A, , k=1,2,...,Nt1 are some time 
independent coefficients). 
Let us look in detail at the case of coupling to a 


Single phonon mode whose frequency we denote Oy: Lneethics 


case the poles of 


& (x) = a (3 29.) 
gra dey! Yi ees ES 
2S 


are shifted to positions 


A+ a ye 
te V¥ (4) +49 (3.10) 


BS 2 


Wemclncecthatetne= probability N, (€) that the atom is excited 


Voemcime t LS egi ven) by 


which is sketched in Fig. 3.1 for various values of Wy: 


N,(t) 


Daleis Sicik “abeye’ jeheleleveleyal lise N, (t) that a two-level atom, cou- 
pled with strength g to a single phonon mode of 
frequency Oye iseexci tedvat timeat;. for, various 


values of Wy- 


58 


= Ve? | ae 


i?.t) 


i. 

aneif.eog o4 Sastif#e szey 
7 "Fea f mA r : 
eo Pe > 2 ee * oe el 
~ 


<*> . (Sos i Afi af” [3) " yi Lisdados; ars 9343 hn28 
“ud cevlip ef 3 gabe 


: 
ae WUE: 5 a 
+ it). “ahs oe = 
a - a ie e 
< y * Bee) ( a 


4 ( =, 


~ >< z “F 
f aveulrey sof I.€ .g69 pt Oaneene al” do Lee 


(238 : 
| - ~ a oe Ps 
{J 7 ae . * i f \, , 

; . es te | hil 
are, \ ; l =* . 

: Lt / 1 /\ 

' Vii AOS 
\ ZZ , f YS 


We note several features: 
ay The time evolution of the excited atom is periodic 
and similar to the behavior of a pair of classical coupled 


; Scharf--, and Mallory" 


pendulums. Recall that Susie 
also found similar periodic time evolution in the case of 
coupling to a single mode. 

(25) ACCOEOM GmeOnud. 93. LO) 7 AG W, * NecleCEcOUDIang 
splits the degenerate frequencies of the two oscillators 
by an amount +g. 


(3) As shown in Fig. 3.1 only modes within a distance 


[Oe = Gg ,0OL 4 gmteract strongly with the atom. 


(4) The atom 'decays' in time ee 1/g. Thus we have an 
‘uncertainty’ relation Aw:T, et hes 
(Se Expanding (3.11) in a power series for small times 
we find 

WhLEs<lpsRg et es G18) 


WoC’ Seg UStmLNeG sLesultm (2.4/7) sof first»order perturbation 
theory. 

With this experience let us now consider the case of 
coupling of the two-level atom to a finite but large number 
of phonon modes. This case can be solved exactly if we 
assume a perfect one-dimensional N lattice-site crystal 
Wien Constant coupling, J. =9 Onsale lien the Nephonon 


modes have frequencies given by w, = c|k| where k=2rn/Na 


k 


is the wavenumber, a is the interatomic spacing and n is an 


integer in the range -N/2 <n < N/2. Since WO, =W_, we have 
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.2atue. Boney 30! 
Y 


N/2 distinct evenly spaced phonon energy levels in the 
interval 0 < Oy < 


PocCgiucnOV PE LNuSHIietiLS Case, =Wwithenompnonons initially 


Dd! where Wy =eT/a is the Debye cutoff 


present, the amplitude 3 (x) is given by 


at SG) a 48 (9) ‘ (3 313) 
N 
x -— A + ne, Ty2 '*) 
D 
where 

Nee l 

N/2 = oe Nx 

205 


The series Ty 2 ‘*) can be represented in terms of the 
digamma or Veeco rtor>, which satisfies the following 
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Using these formulas we can write the series Ty '*) as 


mec = ES = BES 
Tyj2'*) =v tl 20, ee or 
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The approximation in (3.16) becomes exact as the number of 


phonon modes N+. From the discussion after (3.7) we know 


that for any finite number N of phonon modes the poles of 


3 (x) are real. We find the poles Graphically byetooking 
EOtetiicupDOINCSmOEmintcersectzon Of the functions f(x) =x-A 
=ph fe gy) = EY = 

2 Wn N/2 


Notice that the N/2 distinct phonon frequencies and the, one 


(biwithexsreale as Shown ineFiq.,3.2. 


atom frequency have been shifted in such a way that N/2-1 


“yee 


positions that interlace the 'old' uncoupled system phonon 


of the poles (denoted Wir Waren, ) occupy interstitial 


poles and the two remaining endpoint poles (denoted War 


72) have been shifted out of the phonon spectrum 0 <x<w 


D° 
The poles to the left/right of frequency A have been 
shifted to the left/right, with the poles near A being 
shifted the most, indicating maximal coupling to those 
modes. The positions of the poles can be expressed as 
> eae 2n 2 = N_ 
Oe a age + a Oe Tas e Zipeenst ty 5 1 (ere Li) 
where the amount of shift oe of the get pole is given by 
2n 
2 I, See 
NG g N N 
a DN IN si OD Mal! on ) 
Gm COU et (2727 13)) 
n 2 
g N 
\ Ce 


as. Can be sseen, inehld..i> -cerLom scenemLactatiatetiive ee 


cotangent curve is raised/lowered towards the right/left 


by the amount 
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Now that we have found the positions of the poles we can 
invert the Fourier transform using the following integra- 


CHONMCONLOUL 


On N/2 
N wet 
< oa de 2 
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using the fact that 


ace pee a = : ; G32) 


we can write S (t) as 


5 , 
7 [7-2 ay Dy ~1H,t 
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Notercemcne following teatures of (3722)and (3.23). 

(i) iMmeene,expOnentweOlm(o..22)) the shirt OF of the poles 
for large N is negligible and has been dropped because of 
the way the 'new' poles interlace the 'old' poles (see Fig. 
Br )i. 

(2) The coefficient L generally is approximately a 


Lorentzian of width 2ng°N/w, centered about frequency 


C) n aaa 

It is not a perfect Lorentzian because of the logarithm 
term but if the coupling g is small then this term should 
NOEMDE ALOOeimponcaniw. 

(3) Gonverse ly we Gherlast. CWOmCerMS soles Oi22))-, 41s, Owe gLhe 


contributions to S (t) from the outlying poles Wy 


will shortly be shown to be important only when fo is not 
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approximately a Lorentzian. 


heoteusmiebiuserate *hq. @(3e22 )atwith an example where 


La is a Lorentzian and the outlying poles give a negligible 


contribution, namely suppose that 


2 
Pd See eens bet) 8" Spy 
n 2nwy 2) rn 24 2 
(4-7 = At) + (>) 


In Fig. 3.3 we have sketched x, and the corresponding pro- 


bability N, (t) = |s7(t)|i>|? that the atom is excited at 
time t. 
P 
i! 
| 
A. \ 
a Ny Ct) 
! N 
\ _ pias 
I \ An Sits) ih 
! lap hee w 
| AQ= i. 
I 
ge B cls 
p 8 : 
/ \ ce 
D Q 
a 
JS RPEGL eC ae 0 t 
W a 
D Fj 7 TN SG 
= —— ° T= TT =—=— 
DONE ade Real “Tag Long N Dee 
NA' N 
ila 8 eg as, Pa 5 9 temas oA icaeeg 
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1G). mornoren a)e The lineshape f,1 ana (DML ne= corresponding proba— 


Die ey, N, (t) thatethesratomeismeexcited ate time tt. 
Notices that N, (t) is periodic with period Ney yee 
where L is the length of the crystal and c is the 


speed of sound. 
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Notice that N, (t) decays in a time T = 0 / (27g°N) , which 


il 
is just the inverse of the halfwidth Aw of the lineshape 
dn? Thus Ty and Aw obey the uncertainty relation T Ao =1. 


Notice also that N, (t) is periodic with period L/c, where 
Eeisecne) si zeovotmtnencrystal, and) ceis the speed!’of sound. 
This period is just the time that it takes a phonon to 
traverse the crystal and re-excite the atom. Mathemati- 
Calivetices, stemeicmperiodic: because: S1(t) (Eq. (3.22)) is 
given by a discrete Fourier series rather than a continuous 
Fourier transform. (Notice that the perfect periodicity of 
N, (t) is due to the approximation that the poles are un- 
Saurtedmandmunltorm.y spaced. | winereality due to the shifts 
o, we expect Nt) to be ‘almost periodic'.) In the next 


section we show that as the size of the crystal L>+o the 


atom simply decays and is never re-excited. 


Boe DeCayeOrediwAcomeanran [Infinite Crystal 

To solve this problem we can either simply let the 
number, of phonon modes Ngo to infinity in the finite crys- 
talegresihtie( 3622) Or consider the crystal to wbe of anfinite 
extent from the outset. 

In the former case the series in (3.22) can then be 
replaced by an integral, so that the amplitude She) 


becomes 
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now defines the strength of the coupling. 
The other route is to replace (3.13) and (3.14) by 


the amplitude 


ie S 48 (0) __ (3.26) 
ye SUN ae as © de 
or w-x 
0 


where the phonon modes are assumed from the outset to be 


continuous so that 


ete) = oy ee ee = £0) : (6027) 
ie 2 dw 
& Pt eae 
277 Ob < 
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The integral in the denominator of (3.27) results in 


branch points at x=0 and x=wW. which we can connect by a 


D 


branch cut. This branch cut evidently replaces the dis-— 


crete set of phonon poles O47, 050, SO) Of the (finite ’case. 


It is easy to show that the two outlying poles remain. 
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Thus to,sintegrate (3.27) we may use the contour 


(ero) 


On Cy we let x*>x+ie and on Cy we let x+x-ie where 


+ 
See Uenandex sus) nNoweascveal mumber.., Then using the fact that 


lim = = P 
e--0 


[pe TOL) fe (e291) 


where ° represents the principal value, we find that 
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Combining the contributions from G1 and €, into a single 
integral and evaluating 63 and ©, LOueLCmOUELY too moO Lec 
as in the finite case we get just the result (3.24). 
EQes(3224)7 41Smthie exact amp! meudesthatweanwtn heb lly, 
excited two-level atom is still excited at time t, when the 


atom 1s coupledvarbutranrilyestronglyeburceequa lly =to ald 


phonon modes of a perfect infinite monatomic one-dimensional 


lattice obeying the dispersion relation ty IE) NE (eel aus 


study (3.24) ein detail: 


Spee 
essai od euii 


The lineshape 


e(®) = Pe ye REM ae 03.534) 


r 
(wo - A + aan In Fi a (5) 


is to a first approximation a Lorentzian of width [ cen- 
tered about frequency w=A. This is a good approximation 
as long as the lineshape is narrow and well within the 


phonon band 0 <w< w This is clear because the factor 


D° 
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ay 
which makes £ (w) non-Lorentzian becomes large only when 

wo ~ 0 or w x wow. If the lineshape has already fallen off 
at these frequencies then this term is negligible and Z(w) 
is approximately Lorentzian. Fig. 3.4a shows an example 
where <£(w) is almost Lorentzian and hence leads to almost 
exponential decay. If the coupling [ is large, then the 
lineshape is broad and the outlying modes a, and On? also 
become excited to some extent. An example of this is shown 
Timed oms. 4 eee LOmUunderstandsthestime evolution inthis 
case we must use some results from Fourier integral theory: 
If we Fourier transform some given lineshape L(w), then 

the interior of that lineshape gives contributions which 


decay in time faster than any power of t. The endpoints 


oS 3 of the lineshape give contributions that decay like 
ae —iw ts 
se £10, a ye ce" =67 ee And simple poles, of course, give 


non-decaying contributions. 
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Fig. 3.4 The lineshapes and corresponding decays of an exci- 
ted two-level atom under various conditions: (a) T= 
-06 a), A= .5)i lineshape is narrow and well within 
phonon band 0 < w < Op: (b) T= .64,, A= .50,3 line- 
shape is so wide that both outlying modes W6 and 
On /2! also become excited. (c) P= .060,, A=1.050,; 
The two-level atom frequency A lies above the 
phonon band so only On/2 mode and upper edge of the 
phonon band become excited. (The number above the 


outlying pole is the amplitude of that mode.) 


Te 


Thus in Fig. 3.4b there are three types of decay: (1) an 
exponential type decay due to the interior of the line- 
shape, (2) a slower power-law decay with oscillations due 


J . : ' 
to the frequency Wo beating against Oy and Ons2' and 


frequencies oe and w=0 beating against w anda (3) ebeat— 


De 
ing of the outlying modes WG and ny /2 against each other 


with no decay. 
A third possible situation occurs when A lies out- 


Side the phonon band 0<w<w Then the atom is unlikely 


De 
Loedecayer fAthescoupiing Hersenotetoosstrong., An,example 


is shown in Fig. 3.4c where A = 1.05 a, is just above the 


phonon band. We see that the isolated mode On 2 isphaghiy 


excited and the upper edge of the phonon band is only 
slagntily excited’ ggthencontributionsto N, (€) from the out- 
lying pole gives the probability that the atom does not 


decayemm Dewstighnt, oscillation in N, (t) is due to the beat-— 
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ing of the mode w. (which decays like 1/t) against the 


D 
io ole = 


isolated mode w, at beat frequency w 


N/2 beat 


wp/25. 

Let us discuss the signifiance of the outlying poles 
at Wo and On /2° We have seen that these poles allow for 
the possibility that the atom does not decay when the atom 


frequency A lies outside the phonon band. To be specific 


assume that A> ®p- Thensfrompbugs, 342 wegsees that M6 x 0, 


and On /2 is given by the solution of the equation 
Ww. — A-® 
' a ey Ls N/2 Das a i D 
72-4 OF Ltt ora aeas man \ a7 1 ia era : Go -2e)) 


Thus the mode N72 can be identified with the uncoupled 


atom mode A but shifted slightly. Because the frequency 
On /2 lies outside the phonon band, its corresponding wave- 


number iS imaginary, so On /2 represents a mode localized 

in space. In field theory language the two-level atom 
SUbLOUNISE LESC Li awith amlocal cloudsof yvinrtual phonons 
which shift slightly the atom's energy A. Because On /2 is 
a Single isolated mode outside the phonon band, it does not 
decay. 


When A lies inside the phonon band 0 < w < w then 


Dd! 
as we see in Fig. 3.2, the two outlying poles move very 
close to the endpoints of the phonon band and lose their 
Significance as non-decaying modes. Eq. (3.24) also shows 
that their amplitude then becomes small. A band of phonon 


modes near frequency A' becomes excited with a Lorentzian- 


like amplitude 


L008 | a ee reer 
CG) RNar) ech: (=) 
where 
BU ee th 
' Ig D 
eee ue 
A A = ln A ; (237283 .4) 


(Here in the argument of the logarithm we have replaced w-A.) 
and the energy of the excited atom goes irreversibly and 
completely into the infinite number of degrees of freedom 


of the phonon continuum. 


2 


ve 


Three Dimensional Crystal 

For completeness let us study the case of a general 
BSolmencitonalecroyscta li. Althoughetie analysisSutor a 
finite 3-D crystal cannot be done analytically, the analy- 
Shom Oleanmlatinitemo—Decrystal = rollows: that for an 
infinite 1-D crystal with the only change being that the 
phonon density of states is different. Thus the Fourier 


ELanstoumaor Set), namely 


4 ~(x) = eh 5 ; (A, OR) 
baker 
eaee: 
becomes 
g(x) = ae O) ; (30386) 
D 
x - A t+ g°0 pial ee 
0 
where we have replaced 
i 7D 
= \ — | D(w) dw ; (re) 
i: 0 


D(w) dw is the (for now unspecified) phonon density of 
states per unit volume, and 2 is the volume of the crystal. 
As in the 1-D case oS ~ (x) has branch points at x=0 and 

X= Ws (We ignore any isolated outlying poles which it may 
have under the assumption that A lies well inside the 


phonon band.) Thus we may Fourier transform oe (x) using 


the contour shown: 


‘fii 
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| Co eee . 
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complex 


x plane (3.38) 


Along this contour we decompose the integrand in the deno- 
Miia LOreOLm Ss. 50 )mintOuprincipal, and) Imaginary, parts 
AeCOnC tm cORsC an omer. = COMDInIng the Contributions from 


~, and ¢, into a single integral we get 


1 2 
W 
ei ~ixt - < 
Sy (ie) -- | dxe {-§ (x+tie) +d (x-ie)}, (3.39) 
0 
where 
a ee — (3.40) 
soe) a. T (x) 
Waist 2 
7D) 
' 2 dw D(w) 2 
Lae (cet EAR ge Se | tet ’ (341) 
0 
and 
Heys pape (3.42) 


Since only a narrow range of frequencies near x=A (which 
we assume is inside the phonon band) is important in (3.39) 
we may to a first approximation replace Age GOES he CK) = ew 
and ['(x) 2 h(A) =f and extend the limits of integration from 


=o £0 ©. Thus we get 


(ee) i ii 
—-1xt re -~—t 
o~ pa a \entow Vue 
S (t) =5- Sia On bmw wf AAS 


=OO 
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SiO) Gelsiche 
ce - - 
N,(t) = <ils™(t)S (t)|i> = N,(0) e in (3.44) 
i.e. the excited atom decays at rate [f = 27g°D(A)2 which is 


Eiomo cd omit om Gry eeCOUNCED VEL rSteorder. perturbation 


theory. 


Pome NOL cmOUSCNemGcnctal (heony Of Decay, Of Unstable 
Quantum Systems 


Mia 


We wish to present a short summary of the results that 


we, tound in $33.1 and 3.2 on the decay of an, excited ion ‘in 
ENemCONLeXteOfmENe ENCOLY Of G2SSipacion in quantum mech— 
anics. At some time t< 0 we prepare the electron-phonon 
system in a stationary eigenstate of the uncoupled Hamil- 


tonian Ho= namely with an ion in the upper 


dein” “tern 
level and no phonons present. At time t =0 we switch on 
the electron-phonon interaction (Hamiltonian cae Since 
the initial state is no longer an eigenstate of the total 
Hamiltonian A= Be + i the system evolves according to the 
total Hamiltonian to some new state. We are interested in 
the probability that the system remains in the initial 
state, i.e. that the atom is still excited. 

We found that there were many analogies with the 
classical theoryseorscoup led sharmonitcHoscullators..) For 
example in the case of an initially excited two-level atom 
coupled to a finite number of discrete phonon modes, we 


found that the energy of the decaying atom is transferred 


to the finite number of phonon modes and, in a finite time, 


is transferred back to the atom. “This time’ is’ just the 
CIMmen yetakese Lom aephonon) tomtransitethelcrystal. 
Phenomena similar to this have been observed in spin flip 
expenimentecc. In the case of a continuum of phonon modes 
(infinite crystal) the excited atom simply decays and is 
never re-excited because in this case the energy is 


transferred to an infinite number of degrees of freedom. 


The decay of the two-level atom is for the most part 


exponential (assuming that A is well inside the phonon con- 


tinuum). This is indicative of a system that does not 
depend on its previous history, i.e. with no memory. How- 
ever, via perturbation theory, we found the result (for 
finite and infinite crystals) that the excited atom decays 


initially like 
. 22 
N, (t) SC iam 2°) 


1.@€. with zero slope. This is a very general result” in 


quantum mechanics and assumes only that the energy in the 
Bimetoleeccad tem SG mittiiace. me Onethe other hand for very large 
times N, (t) must always exhibit a power law decay «7 1/t 
due to the sharp cutoff of the Lorentzian by the finite 
extent of the phonon band as explained after Eq. (3.31). 
In the rest of the thesis we will ignore these dev- 
iations from exponential decay and assume that S (t) can 


be approximated by (3.43). 
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We would like to point out that our theory and 
results are very similar to those of Weisskopf and Wigner 
on the theory of spontaneous decay. For atomic and nuclear 
systems Razavy and Henley?” have similarly obtained exact 
solutions for the y decay of spinless bound particles. 
Also, Senitzky > has studied various models for the decay 


of harmonic oscillators or two-level atoms due to gener- 


alized loss mechanisms. 


CHAPTER 4 


DERIVATION OF PHENOMENOLOGICAL RATE EQUATIONS 


FOR THE PHONON BOTTLENECK 


4.1 Derivation 

Now that we have studied the interaction of the 
phonon field with a single two-level atom, let us proceed 
to the general case of the interaction of the field with 
many atoms. We begin by rewriting the linearized, coupled 


equations of motion (2.64) for the phonons and atoms, 


namely 
; pike te 
Boi = a + 1 %Szp1e . (4.1a) 
n= = = Sik erate 
== + ~~ . e 
Bor AS L gb, e (4eitb) 


To recapitulate b, is the annihilation operator for a 


k 


~ 


phonon of momentum k and S. is the lowering operator for 
the atom at site r. We are interested in evaluating quan- 
tities such as the probability that the roe atom is excited 


at time t, namely 


t 2) Su ateats See | 
te = Soe Sc 7 (Ane) 


~ 


and the number of phonons in mode k at time t, namely 


ae) = Shs! leisy ee) [Sis (4.3) 


~ ~ ~ 


|i> denotes the initial state of the system at time t=0. 


To solve Eqs. (4.1) we again use the method of Fourier 
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transforms. We define the Fourier transforms of by (t) and 


~ 


S_.(t) as 
hoe = Fb, (t)] ’ 
(4.4) 
Jet) = F[S(t)] ’ 
respectively (see Eqs. (3.2,3) for details), and Fourier 
BoanSLOorme hose 4 i) CO. get 
ne ACO) -ik.r' 
B, (x) = s a J ) e (wees ? (4.5a) 
k era ee xf ve 
i is, (0) eG 
Jj 08) =e wee oy Ene (4.5b) 
Be k is 


We have dropped the subscript k on the coupling Sy assuming 
that the coupling is approximately constant over the range 
of frequencies of importance near resonance. We eliminate 


By. (x) EromeEgs. %4).5b) Sor 3 (x) by breaking the summation 


Over r Sines. 5a), tn toe twoeparts, Anamely 


Seek ae -ik-r _ =e” 
ee eee) E(x) te $ (x), (4.6) 
1g ~ ~ vr'#r ze 
Simo i1 Gala, Jayme COM U4 OD) emeelnuss LG. (4 5b) ais 


replaced by 


E iS, (0) ik-r 1b, (0) 
Ps (x) I e+ Jee Se Oe ) e ae ats = Se ee 
ie 2 z kK X— Wy 
: x-A -) = x-A-) J ~ ~ 
k k k? a k 
5 Lk. (r—-r) 
g € a 
+ ae 2h mace Spies) ete ae 7) 
al x-A-) J ay = “1 
k! C= (i) 


Let us discuss this equation. Notice that the first term 


on the right hand side of (4.7) depends only on the yh 


atom and the third term depends on all atoms but the oh 
atom. From our studies of Chapter 3 we know that the 
fvrsteterm=on) the might®hand side: of (4°°7)" describes® the 
spontaneous decay of the nth AacOMmmeule Leese anitrally 
exe 1ted)-ayoimilaniy, thessecond™term onthe’ right hand 
Side depends on the initially present phonons, and it 
evidentiyadescribess thetexcitation of the (unexcited) re 
atom by these phonons. Finally, by iterating (4.7) we see 
that the third term describes the indirect excitation of 
atom r by phonons coming from other atoms r' (which either 
were themselves originally excited at time t=0 or became 
excited by yet other atoms). The implicit ‘integral' 
equation (4.7) is of course very difficult to solve exactly 
except in special cases. We shall consider these cases in 
€navtensa6. 

For now we Simply derive from Eqs. (4.5a) and (4.7) 
a set of rate equations describing the number of excited 
atoms and the number of phonons present in the coupled 
system. To do this we imagine that at some instant of 
time, call it t=0, the system is in some eigenstate Le 
of the uncoupled Hamiltonian, with several specified 
excited atoms and phonons present. We solve Eqs. (4.5a) 
and (4.7) and find the state of the system a short time 


hom lacer.)  silcesthe time t —.0 LS arbitrary, we can now 


repeat the argument with t=At. Continuing in a stepwise 


manner we can find the complete time evolution of the 
system. 
From our work of Chapter 3, we actually know how to 


find eehestime’-evolution of tthe ren 


atom due to the first 
term on the right hand side of (4.7) (the spontaneous decay 
term) for ta Dletimes;. not justsaeshorte time Ats .Wershall 

see shortly that this is also true for the second (exci- 
tation by phonons) term. The trouble comes from the third 
term because it makes Eq. (4.7) non-local in space and time: 


ia 


It says that what happens to the r atom at time t depends 


on what happened to the other atoms at points ie of the 
system at all previous times. However if we choose the 
time-step At for the rate equations shorter than the time 
tee r—r'|/cethateit takes for a phonon to travel from atom 


ia 


~ 


to atom r then this term cannot contribute to § (x) at 
time At and we may drop it. Thus we may replace (4.7) by 
the equation 

Se OC) WETS iby (0) 


igvig = -—awEa ee Be Ee ee AY 
x-A' + k ik 


~ 


In Eq. (4.8) we have used the simpler expression (3.40) 


Wiehe Avex) co \eandmie> lator EnecmmpLOpaga to Ls 8 (x) of 


~ 


a Single atom in a phonon field. (T, = 1/T is the life- 
HiMmewoOLtlaescing LevotLome muro .64. >) iseiow say local equation. 
LE says that _thevprobability amplitude for the ee atom 


th 
being excited depends only on whether or not the r atom 
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Vere 
a rasaya 


? 


ll aetee Say : 


was initially excited and on what phonons were initially 
presen tencadtetnateacom. Thus: bi cannot describe the 
diffusion of excitations. (We will remedy this shortcom- 
ing presently.) In keeping with this restriction we 
assume for now that the system is spatially uniform. 

We proceed and derive first the rate equation for 
thesphonons. sme liggingnrq 4.8) for § | (x) incomnan (4454) 


EOr By (x) and Fourier transforming we find that 


~ ie 
. = ' a 
Bias ec ee fee 
Demtreore Moab )etig ) elacme SP (0) : 
k k is r a) u = 
rs k 2T, 
i(w-A')t <5 
5) ie a e At , | 
AGE \ 
=a) Det) oo i (4.9) 
Lo a, a0 
a (wow -A'+=—— Cie se Ye ee 
2T, K 2T, 
Let us define the following quantities: 
Q = volume of the optically excited region contain- 


ing the two-level atoms, 
=g is the total density per em? OfStWO- Level 
atoms, ¥ 

N, (t) )=s }<ilsy(t ee Sa tee the densicy ser em oF 
excited atoms, ‘ 

D(w) dw is the density of phonon modes per cm? in 
frequency interval dw (recall that é J+[Dlo,) dw). We 
shall henceforth assume that D(w) * D(A") for all frequen- 


cies near resonance. 


- Tr a 
, 
of fe 
oe = ga i? 
if: ~~ 
L2o0 
7 


apo ty 
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(04) day = <iby (t)b (t)|i> D(w,) da, is the 


k k 


Noh °K 
density of phonons per cm? in frequency interval doy» 


~ 


2 
12 py) Gli = EE ae ala emt Nem hOorentzlan 
(wa -A ENT. 


emission spectrum of an atom, normalized to unity, i.e. 


| £(0) da aris 


Using these definitions we can derive from Eq. (4.9) 
the rate of change of the density of phonons per cm? in 


frequency interval dw at time t, namely 


es f n* 
aE tow Orislel) et (0 Gd) Ny (O)T, (w,t) eS BOS aa pe 
x T, (w,t) (Arsl0)) 
where 
5 -t/T, -t/2T 
Py (w,t) aac (lt+e -2e cCosi(a—A. ye) (Ateias) 
and 
: i 
d t (Sine ee a BER ei; 
T', (wrt) meat 7a area ACOA VAS \ec. Cale 
@—-A' - ~—— 
2T, 
G4 212) 


are rates having the units Secun. ThewrreStetermeonmetine 
right hand side of (4.10) represents the production of a 
spectrum of phonons due to the decay of initially excited 
atoms. The second term represents the resonance absorption 
of the phonons initially present. 

Note that the emission and absorption rates of the 


phonons are actually time- as well as frequency-dependent. 


tog AG mrt as 


a 
4 


siee ugrl't 7 ~amose 
2 ul 7arey 


For example at time tx~0O the phonon production rate due 


to decaying atoms is 


op (rt) = Ny (0) Mee WAY (4.13) 


acy 
dt O 
which is frequency independent because the frequency 
dependence of Ty (w,t) exactly cancels that of the Lorent- 
Zian f(w). Thus phonons of all frequencies are produced 
at the Same rate for time t = 0° in accordance with 
Heisenberg's uncertainty relation AE-At > fi. Notice also 
that the rate actually vanishes at t=0. 

For finite times t>0 Py (w,t) is a rapidly varying 
function of frequency. Experimentally this variation is 
on too fine a scale to be observable, and one sees only 


Py (w,t) averaged over the entire emission spectrum. If we 


perform this average we find that 


f£(o)T, (w,t) do -t/T 
<T, (w,t)> ee on UC Se < Gee 1) 
f£(w) dw 
-t/T 
= — e - < (4.14) 
ik 


For small times this gives the observed rate 1/T, LO 


spontaneous decay. For larger times <I, (w,t)> vanishes 


I 
simply because the supply of excited atoms is used up as 
the decay proceeds. To avoid the latter, we replace in the 


rate equation 
Ae LS 
Ps (w,t) me WAR) ( ) 
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Similar arguments show that this holds also for T,(w,t), 
namely 


<T, (@,t)> => 1/T, (4.16) 


LOGeCIMes. tee Ti - Thus we find that over the coarse-grained 


time scale of order T, the rate of change of the phonon 


density at time t=0 is given by 


NAGE=—0) * Tiere CaN 0) |} 
d =i e a3 oe Ne Ley geen ae 
a Dereon mee (edie) CLO) Ea E(w) do Fa ne S 


(4.17) 


To derive the rate equation for the density of exci- 
ted atoms we simply use the fact that the system Hamilton- 
lan conserves the number of excitations 


) SCE) Same) +) by (t)b, (t) = constant ; (4.18) 
G k 


Or dirferentrating, 


N, + = lil) Gays, (Ee) (4.19) 


vay, fo», Ba 
dtumw2n dts | 


Wrthetgue (4.17 )sethis, implies that 


N. (t=0) 


Sy 


* 


N xs 
D(A") T, pene , t=0) ° (4.20) 


d = =—_— 
—— N, (t=0) = 


at te [au (£() 


To this point we have tacitly assumed spatial 
uniformity. The phonons were non-local plane waves 


extending over the crystal and we replaced a }<s7S "> by 
ie ~~ 


the average density N Let us now assume that phonons 


3° 


are absorbed or emitted locally at the position of the 
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atom where the interaction takes place and allow a small 
amount of spatial non-uniformity in the quantities an and 
N,- PpeneECS a (4a) ands (4.200 )aresstidl ‘locally valid. 
With the following modification they are valid everywhere: 
We add a phenomenological diffusion term Bo) V" ny (eyo, t) 
to the right hand side of the phonon equation because the 
phonons tend to diffuse from regions of high to regions of 


low phonon density. For the frequency dependent diffusion 


constant %(w) we choose the classical expression 


ua) ae 3 = yay (Arb) 


where A(w) is the mean-free-path, 


il N* 


i (ea D(A) cT, 


£ (w) (4.22) 


is the absorption coefficient per cm of phonons of ancular 
frequency w, and c is the average speed of sound. 
If we now assume that these equations are valid for 


all times we get the equations in the final form 


ie DSH saarOd se eekanida i Sey wen eae a(w)c 

2 

' + Blo) Von (rio, t) do , 

| 5 N, (r,t) 

lie N, (r,t) “1.5 See + | au Ma geo s Gua). Com . (AG23) 
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4°25) Comparrson with Pauli"s Derivation of the Master 
Equation 

It is instructive to compare the assumptions that 
went into our derivation of the rate equations (4.17) and 
(4.20) with those that went into the derivation of Pauli's 
master equation or Boltzmann's equation) -'>~, 

In all cases an important assumption is the existence 
of distinguishable time scales, an idea first introduced 
by Bogoliubov>>. The first characteristic time Ty is the 
Spode GeO mmcniecwiuraclon of achesinteraction) of a phonon with 
a Cr? ion. For times t<T, we saw that in our rate equ- 
tions the rates Py (w,t) and PT, (w,t) change very rapidly 


and a detailed quantum description of the system is needed. 


ioeciesekinetrceregqume, i.e... 10r times t> T- the rapid 


1 
variations of Pr, (wt) and Py (w,t) can be averaged out 
leading to the kinetic Eqs. (4.17) and (4.20). The second 


characteristic time is the mean-free-time 
T, = A(w_)/ce = ——— (472 4) 


between interactions of the resonant phonons with Cron ions. 
Roretimes | ty T. the resonant phonons will have interacted 
with the Cru ions several times and reached a local equili- 
brium. Thus in this hydrodynamic regime Eqs. (4.17) and 
(4520) sare slocabiyavalid mande)teusmsureicirent: COsaddra 


hydrodynamic diffusion term to get a set of equations (4.23) 


which are globally valid. 
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In the case of the Boltzmann equation, the assump- 
CLoOnetiia t Ty ce T., MeauS ethatelten smrestricted to, the 
deScu pElOnmOLlauciiute gas... ineoumgecase, using, Eq. (4524), 
Wemcccmtiiad bart amcancechatecthe deisitymo f Cran ions in the 


E state must be less than the density of phonon modes on 


resonance, i1.e.: 


N*<p(A) Aw. (4.25) 


where Aw =1/T, is the Halfwidth of the resonance. This 
restricts us to rather low densities f 2 ie eee However 
Since the Boltzmann equation works well at much higher den- 
Sities than one would suspect, we hope the same holds for 
Our rate equations. 

To derive the Boltzmann equation a second assumption 
is required, namely the Stosszahlansatz, i.e. the assump- 
tion that each collision destroys the correlations that 
have built up between the two colliding particles during 
the time since the last collision. This allows a dynamical 
description of the system in terms of the one-particle 
ic ieeroe Onmrunct1on fy alone, which introduces irrever- 
Sibility. In our case and in the case of Pauli‘s master 
equation the Stosszahlansatz allows a description of the 
system in steps of time At where T,<At<T,- In each time 
step the system begins in an eigenstate of the uncoupled 


Hamiltonian without any correlations (e.g. in our system 


only the number of excited atoms and phonons is specified. 


This 1s a frequent criticism of Pauli’s derivation of the 
master equation.) In the case of the master equation, the 
eho. ce T, <4t<T, allows one to find the transition rates 
UsiIngmelustweovder pemturbation theory. | Injour case at 
UuULOWwWSmMUSELOSGLTOD ethesnon—lOcal termmein u(428/) 2 In both 
cases the result is a set of first order differential 
equations which describe the irreversible evolution of the 
SySctem': 

These parallels between the master equation and the 
bottleneck equations are summarized in Table 4.1. 

In conclusion we emphasize that the bottleneck equa- 
tions (4.23) are valid for densities N* low enough that the 
CC nC ise On T, < T, (N*) holds. At higher densities where 
T,(N*) < Ty we expect that we must include the non-local 
term of (4.7) in the rate equations, yielding a new effec- 
tive lifetime T, of a local cooperative group of two-level 


i 


atoms. 
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CHAPTER 5 
SOLUTIONS OF THE RATE EQUATIONS 


5.1 The Solution and Comparison with Experiment 

In this chapter we solve the rate equations (4.23) 
describing the time evolution of the coupled system of 
phonons and atomic two-level states of excited eae ions. 
Weebegtn=by rewriting Eqs® (4/23)hin atnotation more “in 


line with previous works on the phonon bottleneck, namely 


3 Nye ee 
rs ple Cote Key = ielepiahy Sema Ips Chios sys ava (vu) ¢ 
+ Biv) Vn, (x,v,t)av , (5. 1a) 
phe = 
loge nen? Poy a 2 [av nog (ery et) atv) : (ee) 
Beh ere) Cy is the density of phonons per em? in frequency 


iotewavaledymaboubsy —w/27 and at position r, and N, (r,t) 

is the density of em ions per om? inethem Aes ta tomar 
position r. The fuestterm= on thewruqhte hand side of s(Ss-ia) 
deceribes the increase in the phonon density at position xe 
and frequency v due to decays of Crom ions in the 2A state 
at position © occurring at the yrate 1/T,- Due® too these 


decays a spectrum f(v) of phonons is produced. We assume 


(for now) that this spectrum is a Lorentzian 
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with full-width-at-half-maximum AVG centered about the 
resonance frequency Vpeteeall andenormarezed=to unity. The 
lifetime Ty of the en ion against phonon emission is T= 
(2mAv.) >. We regard T) as the single adjustable parameter 
of the theory. 

The second term describes the absorption of phonons 
by atoms in the E state. The frequency-dependent absorp- 
tion rate a(v)c per second is peaked at the resonance 


Erequency. Wevassume that the absorption coefficient a(v) 


per cm is proportional to the emission spectrum and given 


Bye? 
a n* ae Ne 1 
ek pCa 4D CLL Ay oe) (5-3) 
1 fe) 1 ef ©) 
Av /2 


Here D(v)dv is the density of phonon modes per em in the 
ruby crystal in frequency interval dv which we assume to be 


the Debye spectrum 


Buy ak = DAE R RSS) eas (5.4) 


127v-av 
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Q 
Boe 


oa 
c? c 
My 


ator Jie 10m cm/sec, n= Wes aon cm/sec and aa oe Mal 


cm/sec are the transverse, longitudinal and average sound 
velocities respectively of ruby??. Our use of the Debye 
spectrum is justified because aes ee The quantity 


s T “ T : 
Ds a) DIC) 7 AVG x D(v,) x Av is a measure of the density 
3 


O 


of phonon modes per cm~ "on speaking terms" with the decay- 


ing atom. (oe is the width of a rectangular spectrum 
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with the same area and height as the Lorentzian (5.2). The 


poo ioaC COs Nestor lelsseri xed by thescond: tion enees” 


N* 
D(v.) Cc T) 


(ary 


We can understand the expression for a(v) as follows: 
According to equipartition of energy a quantum of energy 
Avy added to the atom-phonon system is as likely to excite 
one of the N* atoms as one of the Do phonons that are all 
at approximately the same energy Avo and are all approxi- 
mately equally coupled to the N* atoms. Thus if the life- 
time of an excited atom is T, seconds then the lifetime of 


1 


a phonon of frequency we must be Tx (Do/N*) seconds. 

Since this phonon travels with velocity c, its mean-free- 
path is A(vg) = 1/alv,) =D eT /N". The absorption a(v) per 
cm of a phonon of any frequency v is reduced from a(v.) by 


just the Lorentz factor 


al 
Va 2 


O 
ibe 72? 


Hence we wet Eq. «(5.3).. |For ruby we have Div.) =7.9 x 107 


sec/cm> 

The third term describes the decrease in density of 
DMiOnOnNS*ac  pOSTCION ceo y edi fuSlOneOLmpnonons saway from 
position r and eventually out of the excitation region. 
The frequency dependent diffusion coefficient @(v) is a 


minimum on resonance and increases as |v-vo| increases. 


We take (v) to be given by the standard kinetic theory 


expression 


SA eo Oe ae (5.6) 
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where A(v) is the mean-free-path and c is the average 
phonon velocity. Eqs. (5.la,b) together implicitly take 
into account the fact that the time taken to traverse a 
mean-free-path is the sum of the time between collisions 
and the time T, of the collision itself. (Note that Ov) 
actually increases without bound as Pare oS but this 
pathology will be remedied shortly.) 


Finally we impose the restriction that Nj,(t) << N* 


2 
for all times so that we can neglect reductions in N* due 
to excitations as well as neglect stimulated emission. Also 
we neglect the continuous, slow generation and loss of ions 
in the 2A and E states due to optical pumping and Ry and 


R, Gadlatlona | NOting that the lifetime of a longitudinal 


phonon against anharmonic decay is ~2 us and that of a 
transverse phonon even longer by orders of aeearenste 2°" 
we neglect anharmonic phonon decay. We defer until 85.2 


the possibility that the spectrum f(v) may be replaced by 


an inhomogeneously broadened Voigt profile. 


Spatial Dependence of Solutions 
EoS2esto.la, Dimarewcoupled fine tseorder, Jtinear, homo- 
geneous differential equations and can be solved as follows. 


We group together the phonon source and sink terms as 
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N, (r,t) 
Sve, Upon) = Getty) Ch Sere -n 


Gap FOP ial ewer Oe keer (yeas) 
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ph 


and write Eq. (5.la) as 


) s 2 
aid Nop (ErVet) SS Sage) Ae H(v) Von (x,v,t) dune (5 8) 


ph 


We try a Separation of variables solution of Eq. (5.8) of 


the form 

wage Medic = Non (vet) dv R(r) - (5.9a) 
and 

N, (r,t) = N, (t) R(r) , (5.9b) 
so that 

S(xr,v,t)dv = S(v,t)R(r) 5 (529C) 
didevetic ott and N, must have the same spatial dependence can 


be seen from Eq. (5.1b) and is essentially due to the homo- 
Gener tyaOrehdS (ola) rand (5.1Lb) 
SUS emc Clones milo. 9) wintoulEd.. (5.8) (gives usethe 


two equations 


d 2 

qe ge teehee (ua biiec ple siersntlyares F (5.10a) 
2 2 

Ves RECT). aes (alate Cs LOD) 


where ae is the separation constant. If we assume that the 
excitation region (defined by the pumping laser beam) is 
atelicin bcency lndemweOrsbadlus is=ib/ 2s (and) that we there- 


fore have cylindrical symmetry) then the solutions of 
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¢t lens 
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EOtme os lLUDjeare the zeroth order Bessel functions R(r) = 


Ue (eee) eet ne=bOoundary Condition that N, (and therefore! 


fe) 2 
ety vanish at the edge of the cylinder r=b/2 requires 
that the separation constant k takes on only the values”: 
a Biya TNL at alliance (5.11) 
b y i SOE | tay. aaa : 


Here x is the ay zero of Js) - The general solution of 


Roem orac)) WeLSm Elen 
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co K 
an, Oa Date eh =) CiI (Kp eF) ny iy, t) dv (5en12)) 


ele i=l 


where the c,'s are constants prescribed by the initial 


1s 

conditions and non (Yet) dy is@the solution of Eq. (5.10a) 
K 

with mie The mode oan is the only mode excited to any 


appreciable extent if the initial condition of a spatially 


uniform laser pulse obtains. Furthermore, by inspection of 
K 
Eq. (5.10a) we can see that a is the mode with the slow- 


est decay rate and therefore the one of primary concern to 
us. It is interesting to note that Ky is rather insensi- 


tive to the precise geometry of the excitation region. 


Nor: course there are phonons outside the excitation region 


but they can no longer be scattered back into the region 
to participate in the electron-phonon interaction. There- 


fore we can take Pe AY raala fees leys 7 
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Time Dependence of Solutions 


Let us now obtain the time dependent factors of the 
K 
densities, namely nop (vr £) av (henceforth we omit the Ky 


superscript) and N,(t). They are the solutions of the 


coupled differential equations 


d No (t) 
enon dv i=sti (adv T; ~ Non (Y7 2) av ov yc 
eye ae Oy My Ft) dv (5.13a) 
t4 tv) ph =~ 
d Ne 
aE No (t) =- i = [av Magy es CMG) xe (Seat 5.) 
where 
t.(v) = Genoa. sat (5.14) 
d b si 


is the diffusion time of phonons of frequency v. We have 
noted already that (5.14) must break down eye | en large 
enough that the mean-free-path is longer than b, because 
then (Vv) >> be and tg) D/C mt om se UlaeCepLaDaLe 
because the minimum time required for a phonon of any fre- 


quency to leave the excitation region is the ballistic 


escape time i, Siayel Thus we replace (5.14) by the gen- 


eralized escape time 


? 
z ies al b 
Ce) = ty + tg (v) as + corey Re ((syealisy)) 
Then for $(v) >> bce we have ea~ t =b/e and etOre. (oj) <a DC 
we have oe «ta: SOlwingeigds. 05. la )eeoy = the=method of 


Laplace transforms we find that 


eee | st 

a ae, oi ° dse me Ou ; (Ec ise) 

and © 
peng” k st 
N, (t) mh Spee } dse N, (s) ; (ore Ob) 
4 
where ae aaa and N,(s), given by the Laplace transforms 
Que (Hep sees ae wale 
~ £(v) 
{nv t=0) + rT, N,(s)} 
ls MOD St Se ee oe a a C57 ai) 
ph 1 
SET ONY) Coch os mi 
e 
with 
f = - pee ber (A) Cee 
iN, (t Oth ae Jf Cah Let) State +1/tew)’ 
N,(s) = 
2 e 
ue, oes’ S,) 
Ty (Sil D) 


ee ee a Se 
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‘Thesaverage: Dallisticupatheleng thm. ainsidesa, cylinder got 
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radius am andvlengthie- res approximately <t2=r 4 (1n=$-1). 


r 


Tneour case. =) s35emm and |, = L0amm so that. <2> ~ 205= b. 
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where 


Lis) = aly ely) (5.17c) 


st+ta(v)c + 


| to fv) 


Genera lalecabitness Of SOluEiLons 
The contour © encloses the various poles and branch 


cuts of n v,s) and N.(s) which are located in the half- 


ph 22 
plane Re{s} <0. Let us begin by studying the singulari- 
ties of N, (s) in the simplest case when there are no 
phonons initially present. If we replace the integration 


over v in )(s) in the denominator of N,(s) by a summation 


2 
over the phonon modes, we can locate the singularities 
Graphically by@iinding the valnuesgot s for which = le =e) (s)\- 
As we see in Fig. 5.1, ) (s) has a pole due to each phonon 
mode included in the sum, and the poles lie spaced in the 


range s, < s < -c/b where 


R 


bo 
one 


jie = [avg ‘ - et ee (5.18) 


jh fy Re. Se ered oy 
(1, BP 5 


The pole at s=Ss, is due to the resonance phonons which are 


R 
the quickest to disappear due mainly to resonant absorp- 
tion. The pole at s=-c/b is due to phonons far-from- 
resonance which disappear relatively slowly via ballistic 
flight. Coupling the atom-phonon system causes these poles 
to shift to new positions determined by the roots of the 


equation -T\s = ithe The graphical solution shows that 


all but two of Ehe roots are located in the interval 
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Sp sss -c/b and are shifted only a small amount. One 


s,o7l 
1? , has been 


shifted to the left of -T)" and when inserted into Eq. 


(5.16b) yields*the initial fast decay in time T 


pole, whose position we have labelled s=(-T 


' 

< 
i Ty OL 
the coupled atom-phonon system. The other pole whose posi- 


tion we have labelled ee, has been shifted into the 


interval -c/b<s<0 and represents the final slow decay in 
time T. of the coupled system. The poles in the interval 


S, ss < —c/b, when inserted into Eq. (5.16b) yield the 


intermediate decay. We should point out that Fig. 5.1 as 


drawn describes the situation for densities N* ~ 10+4 19! 


cm for ruby. As noted in the introduction, at lower 
densities the phonons escape ballistically, a(vo)b geo de 


and the interval Sp sss -c/b becomes narrower than Fig. 


5.1 suggests. On the other hand for densities NagenD San 


O 
Ome a=, the resonant phonon lifetime is less than the 


2A state lifetime ue and -T)" 


interval Sp es a —C/ De 


Let us now go back to the continuum limit of } (s) 


actually lies inside the 


and replace the summation over phonon modes by the original 
integral. A detailed calculation shows that the two poles 


located at ms ai and s=-T. remain but the set o£ 


i 


eoles in the internvalmsc™= <)Sec-C/beismrcD laced iby. a 


R 


branch cut connecting branch points at s=s, and s=-c/b. 


R 


Thus the only change in the above description is the 


rather minor change that a discrete set of intermediate 
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decay rates is replaced by a continuous one. Indeed we 
find that the time evolution of the density N, of atoms 
in the 2A state, with no phonons initially present, using 


HC seen LOD i mese dl VEN by. 


Ny(t) | Seay ey is =) -t/T, 

N, (t=0) ~ as uy dso sods) + Re (52193) 
S 
R 


where £(s) is the discontinuity across the branch cut, 


f(s) = lim st { st + — \, (5.20) 
e>0+ + = )(s+ie) s+ )(s-ie) 
T 
i 
and Rie Re are the residues of N, (s)/N, (t=0) at the poles 
3 = siya, Sl eprespectively)§ io. 
R_,R z (5.21) 
te Ss i d Tv 
ne ee ee 
u sat panes 
TS Ts 


These quantities of course obey the normalization condi- 


tion 
-c/b 


Lar se | No(s) dsa=al a (5.22) 
S 
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Solutions of the Rate Equations 


Fig. 5.2 shows a typical example of the time evolu- 


tion of the excited atom density N,(t), given by Eq. (519), 


showing the fast initial, intermediate, and slow final 


decay. Also shown are the total phonon density Se! = 
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N, (t) +N 


(v,t), and the total density of excitations 
al in the region. 

Let us now look at the final slow decay of the 
system in detail. We recall that its decay time is Te? 


-1/s where s lies in the interval -c/b<s<0 and satisfies 


the equation 


(oe) 


-T)s = ) (s) = dv £(v)il-Po(v)} , (5228) 
where 
&: CHEV) Ke: 
Pi (y) claatshaibina:. Ada Eres Sa ee (5.24) 
Sara (i) Cot EGR) 
e 
In general Eq. (5.23) must be solved numerically and these 


results are shown in Fig. 5.3. Here we plot T. vs N*/D, 
for various values of T,c/b. Note that Doe and thus the 


horizontal scale, depends on the time T, of a particular 


1 


curve because of the relation 


_ Pie / 220% ag! Baer 
D — == ae e 
fe) aT, Ty 


Fig. 5.3 clearly shows that there are two density regimes 
with quite different physical processes governing the 
decay of the system in each. As we will see in detail 
below, for low densities the system decays in time deni 
jeywe (al cs a(vo)b), which is the ballistic or free-flight 
time plus -thewresonant phonon diffusion time. Because for 


large a(vo)b (i.e. large probability of phonon scattering) 
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Solid lines: Theoretical final slow decay time Te 


as a function of the density N*/D, for various 
values of T,c/b, fromebqss (5523). Notes that the 
order of lines is reversed at low densities. 
Dashed line: Experimental results from Refs. 20 
and 21 (see also Fig. 1.5) plotted assuming D.=10"7 
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this would suggest that Tye Ib we may call this the 
Spatially diffusive regime. However the figure shows that 
beginning already when a(yo)b ea Ae T. hardly increases as 
we increase N* by, OLCeCrs OlLemagnmcude,;mwindicating a 
process competing with spatial diffusion. In the high 
density regime a(vo)b >> 1 this process, which we call 
spectral diffusion and for which we will find that T x 
T,¥a(vo)b, completely dominates. We also show the decay 


DOnncas 


times found by Pauli et al (Seca bigs |. 5) tor two 


different excitation region diameters b. For these two 
experimental curves Do x ome a was found to give the 
best fit to the theoretical curves so.that T) a MM) Fall 


We shall say more about this rather large lifetime in 85.3. 


Decayelimesr  inecase N4 << D 
tS 


For the limiting case Now Do we can actually find 


T. analytically as follows. We assume that the absorption 
* 
PueoObabwiie Ly ac. () =) bee WW ee iL hele) iegiske Talsl Gbichyleteeue ley fens 
fe) Doct, 
the excitation region is so large that T, eee Vey he © algh falas 


ballistic regime we cannot speak of true diffusion and we 


replace the diffusion time (5.14) by 


tay) = Ses) (Sa P a) 


which is the ballistic time multiplied by the (small) 
probability of the phonon being absorbed and scattered in 
this time. Keeping the first three terms of the binomial 


expansion of the reciprocal escape time 
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-l 
aay = (20 tab), (5426) 


Eq. (5.24) becomes 


Gores) 


endethe gntegqrand of bq. $5.23) cangbe. written’as a sum of 
powers of Lorentzians. These can be integrated exactly so 


that Eq. (5.23) becomes 


SOV. SS ee 3. 94 SS ox Gamo) 


Since the right hand side of Eq. (5.28) varies very rapid- 
ly for s=-c/b (see Fig. 5.1) and the left hand side is 
relatively constant, we replace the left hand side by 


-T,s x T,¢/b anc maicmecheomscOlUtLONnmOLehd 15.20) tobe 


a(v_)b 


(LS See 2 boc) ; (5.29) 
2(1 - Ty 5? 


The final slow decay time the Ou N“/D, << 1 is therefore 


Ca Qe) 
Cha eee 2 ona & (5.30) 
2(1-T, >) 


Thus we get the reasonable result that the system decays 


in a time des Se ABs sta OSs , which is the time an "average" 
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phonon, namely a phonon a half-width from resonance would 


escape in. 


Decay Time T. inmCaseatl a =D 
nn —{ ) 


In this more interesting case we assume that the 
absorption probability a(vo)b >> 1, so we are well within 
the diffusive regime. Since 1/t, (v) << c/b and we demand 


that |s| <c/b, we can approximate (5.24) by 


v-Vv 
Cc fe) 
1 when COV ORS iv, 72 < va (vo)b ; 
Pv) = G55 1) 
a abe oh 
Lo when a(vje <p <> iy,72 | > Ya(vo)b és 


SUpStUucuUcINGm (Deo lnto. bo. (5.23) ewe get 


-1,5=2 | 20 nd yy ee (5032) 
1l+x Ya(v)b 
va(vo)b 
where 
v-v 
x = 2 ’ 
Av /2 
or 
Tig eee array (5633) 
S S oa | re) : 


This is just the decay time found by Hote ceana and 
Veklenko/ by rather more complicated arguments. Eq. (5.33) 
is a very interesting result. It says that the final slow 


decay has absolutely nothing to do with spatial diffusion. 
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Indeed approximation (5.31) requires only that the absorp- 
tion probability aCvojb >> 1 and would have held just as 
well had we chosen the diffusion time ty(v) =0. (Of course 


we must still have the total escape time tv) Pa. Cre) 


In other words when a(vo)b Fels pat ale dL erusvoOne: Ss SO 
Slow and ineffective as an escape mechanism that, as we 
shall see in detail below, the phonons use a shortcut - 
namely spectral diffusion. By this we mean the following 
process. 

Excited atoms in the excitation region decay in time 
T) emitting a Lorentzian spectrum of phonons. Those 
phonons in the wings of the Lorentzian satisfying 
Co De alls (1). Compare ractaonwie= 2/(mVa(v_)b) ) escape 
unhindered. The rest are absorbed in times t << T, and 
Subsequently re-emitted again in a Lorentzian spectrum. 
The process repeats itself again and again until all the 
phonons have escaped "through the wings" of the spectrum. 
From this simple argument Eq. (5.33) follows immediately. 


Although some spatial diffusion obviously takes place 


concurrently, we see that it plays a secondary role. 


Phonon Spectrum 

To prove these assertions let us now look at the 
time evolution of the phonons. We recall that Dea) dy 
1S givensbys (S2iGa )mwithes(5% 1+). ein general these equa- 
tions must be solved numerically and some typical results 


are shown in Fig. 5.4 where we assume that there are no 
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Fig. 5.4. Time evolution of the phonon spectrum as Es 
for case Ny Det and b/(cT,) = 3. Note that as 
t+o the spectrum is approximately rectangular 
with width v= Ya(v yb Av: 
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phonons initially present and that the absorption probabi- 
icy, a(vo)b is large. We see that for very small times 
some of the excited atoms decay and emit a Lorentzian 
spectrum of phonons. Soon the resonant phonons are 
absorbed by other atoms causing self-reversal of the 
spectrum. Since only the off-resonant phonons exist for 
an appreciable amount of time before they too are absorbed, 
the spectrum broadens and flattens. The spectrum continues 
to broaden until time t=b/c. By this time the far-from- 
resonance phonons are escaping ballistically and those 
now-empty modes can be re-occupied by other phonons whose 
excitation energy comes from phonons lost from the middle 
of the spectrum. Thus a steady state is established with 
a phonon spectrum which is continually losing energy from 
the center to the wings and which decays away uniformly 
with a single decay time To: 

This slow decay of the phonon spectrum can be 


deri vedeanalyticallyey, BromyEq. (>4l/7a), we-haves that 


EQ al 

n y f¥ 1s) dv = D(v)dv - I : < (5i4)) 
P sta(v)c + N 

ie thy 

S 

es 
Pier factore (ser onG ica Ee has a pole which lies in 
e 


the interval Sp <s<-c/b (see discussion before Eq. (5.18)) 


and the factor N, (s) of course has poles at s=-1/T, and 


Spt edly fale and a branch cut in the interval s,<s<-c/b. 


R 


The final slow decay therefore is given by the residue of 


La 


Nga Sa Phe aceathe slow=pole"! s = -1/T.- We find that 


N, (t=0) Bey abe 
De eee yet) Ca Di (0k) ce (v) -—;-——_R e 2 
parce DD O igh Aah N SL 


(eo 5)) 


Thus the final broadened phonon spectrum, aityar Cop aly 


does not change shape but decays uniformly in ee Te as 
we asserted. 

Recall that torslarge densities N*~ >> D, the spectral 
shape Saat UseOUeoebVatnenaporoximatlons (5.31) ,.1-e-% 
titers equite flat. andphas. breadth, Av = Ya(vo)b Avj- The 
slight self-reversal is due to the negative term s=-l/T. 
in (5.24) which is relatively important near the 
"Shoulders" of the spectrum. 


For low densities Nae Dy EOS 005.275) ands (a7, 29.) 


give the phonon spectrum 


Cc 
a BD ; 
Peay hae ee Zo oie = eter 
( AP {1 eo Sa Se 
Romy/2 L VET } 


2 
O 
Ube ) 
Av (/2 
whose half-width appears to be increased from the intrin- 


2 


but which is equal to the original Lorentzian in the far 


Sic Lorentzian emission spectrum by a factor x (l+3a(v)b), 


wings. 


Ii3 


Quasi-Equilibrium Populations 

Finally let us look at the quasi-equilibrium popula- 
tions of the two-level states and the phonons. We now 
generalize and assume that the initial laser pulse excites 
both phonons and en fons in the 2A state. Then Eqs: 


(5.16) and (5.17) become for the final slow decay 


Nace) N, (eq) -t/T 
: =r ac 7 (5.37a) 
n* n* 
where 
* a Se ’ CSres 7D) 
N x 
N ao yay NY, | 
and 
ny, (v,t) dv N, (t) 
Diy) ava oF ici Sean ora (5.37c) 


Evaluatinastherrelativer2A state population in the limit 


Nae D, we get 


ee [Vey acpen, Cr 
N, (eq) q N. (t=0) trecaeeae (Th Sh Hentd a(vo)b + esos) oo 
* * * Cume5 Hie ’ i e 
N N + 2D (1 - 2T) HT Doty) 
where 
te SO Oe eee (56s ob) 
fe) O82 sephe so 


is the density of initially injected phonons in the fre- 


1 


quency interval AVG x: 


= bee tt 

at Ou «).2 

x — — = = Se illinois, 

7 aur de 
ae BY f 
a 
“23 i exonot 
> tJ 


< 
and 
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Summing over phonons in (5.37c) we get the total 


phonon density 


Ng lS) = | n 


(OC) 


MiemiactOLeOLmeealneLront.oL nae and De thigy Sig ehsken) yee! 
eee eee erom ncetollowincmract-meTf aecr-1 ion 
absorbs a phonon (which it does according to a Lorentzian 
of distribution), another phonon will be re-emitted with 
certainty. On the other hand if a phonon is emitted with 
eeeC rent tiaied obit butl On ,etnens 1 tacubsequently may be 
absorbed (again with a Lorentzian distribution). Since 
thceprobabilittyeO1s the latter process involves the product 
Omeewowhorenuzians, witheintegral half that of a single 
Lorentzian, the phonon population is favored by a factor 
of 2. (Indeed if we replace the Lorentzians f(v) and a(v) 
by rectangular distributions the factors of 2 do not 
appear.) As one would expect, increasing the lifetime Ty 
or the absorption probability a(vo)b Savoremeneeenectates 
in the quasi-equilibrium. 


In the limit N* >> D. Eq. (5.37b) becomes 


he rt=0)Va(vo)b a 


n* N* +D(v.)¥aW,)b Av, 


N,(eq) N.~(t=0) +n 
a fe (5.392) 


In this case all the initially injected phonons and phonon 


modes in the broad interval Ya(v yb )b Av, (7.e. those 


- 1 
7 
2 ~ > ; » 
w 

i ‘ry | 

. 

The. 
~<, 
- 
‘ 
a i 
x ~ ~~ ' : - ~~ cy 
‘ 7 b 


‘ 
i 
= , ] 
> ‘ . 
~ 
« ad 
> + tei f, 
¢ ‘ “+ -s 
— 
4 a *% 
. ° 
. -_ 
~ " +. J 
— ( 
7 
| P 
~ ’ 
_ = 7, 
do 
- . rr heel % 
. ; 2 ! 4 
= y 
~ ‘ee -¢7; aa 
’ f re a ‘ = | 370 
a 
— . 
On 
* 
¢ 
: i > 
e - F] . sz “ye 
¥ ¥ ' M« sf J £ J oa | F ee 
, 
ree 
‘ - 
¥ 7 bai 4 
é 


: - : - ' wt aor , wig2edsa sna 


_ ¢ 
ay = eew 
os 


26-£ E205 
nll : 7 
gil ois & 


yan < 


Saulo my UngetnescOnd tionvd (Vy) bees) participate in: the 


quasi-equilibrium. The total phonon population is 


Nee) 
=se24 feral (8) 20 Wal 2 ‘ (5693.9'5)) 


WERSOMVCME Cano oD) MmNUMNeEr Ca byaand: Nn. Fig. 5.5 


plot the quasi-equilibrium population N. (eq) /N* vs ae 


2 


for various values of T,¢/b assuming that no phonons are 


initially excited (n_, (v,t=0) =0). We see that for N*<<D. 


ph 
(see Eq. (5.38a)), No (eq) /N* =N, (t=0)/ (2D) which is 


independent of N*. As N* increases, the phonon band in 
quasi-equilibrium with the atoms broadens like Av = 
Va(v Jb Ay, sO that in this regime (see Eq. (5.39a) 
N, (eq) /N* = N, (t=0) /(D va(v_)b ) which varies like N* ?. 
Hina blye LODEN = > Do the N* term in the denominator of 


(5.39a) dominates so that N, (eq) /N* x N, (t=0) /N* which 
varies like Ron we 


These calculations can be compared to experiments of 


Pauli et Aisa ae if we note that in those experiments the 
: Ofte 30 , fi, 
initial laser pulse does excite phonons as well as cr? 


ions. Indeed in their density regime ay 22 eal 


La aie the appropriate curve, namely T,c/b =1/5, behaves 


SN eS 


like 


~4 


N, (eq) /N* c N* (5.40) 


in excellent agreement with their experiments assuming 
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Fig. 5.5. Quasi-equilibrium population N. (eq) /N* vs N*/D 


O 
for various values of T,¢/b assuming no phonons 


are initially excited. 


5.2 Inclusion of Inhomogeneous Broadening 

Until now we have assumed that the absorption and 
emission spectra are described by a Lorentzian profile 
centered about the resonance frequency Us However random 
static strains in the ruby crystal (or higher temperatures) 
can cause the resonance frequency to deviate locally from 
Vo: This results in so-called inhomogeneous broadening. 
The central limit theorem tells us that for large statis- 


tics the deviations have a Gaussian weighting. Thus we 


can replace the Lorentzian by the Voigt profile 


00 2 
Av 4 
f(v) = atl e.9: = 5 ; (5.41) 
Ven) = eam + (1p 1—2y:) 
where 
vrv 
n =2 Vin 2 2 (5.42) 
Av 
G 
and 
AV, 
Be SRE Aben 2 oy (5443) 
Ave 


AVG measures the amount of Gaussian smearing of the ori- 
ginal Lorentzian of width Av: As before f(v) is entered 
about a andenormal’Zedatomun cy mush OOmCONnGltGlOnm Uo.) 


we can write the absorption spectrum as 


We PA * oF 

~ os 

ND) 2 we Gao) | aay “4 (542) 
Oss leeeieey i Co ee” Se (ae) 


Note for Av, = 0 (Av = 0); f(v) is pure Lorentizan 


ey, 


Cass lan) MOSt@mimpontantly, £LoOrefinite AV /AV or f(v) is 
approximately Gaussian near Weis and approximately the 
unbroadened Lorentzian in the wings. 

We can now apply our method to the Voigt profile. 
For, example solving Eqs. (5.23) and (5.24) numerically for 
the slow decay time ite yields the results shown in Fig. 
5.6a. Here we have plotted T. vs Nas (ere T,c/b = 1/200) 
for various values of Av/Avo: In all cases for small 
N*/D. we see that the decay time T. is smaller for the 
broadened than for the unbroadened spectrum simply because 
the broader the spectrum the smaller is the absorption 
a(vi)- For somewhat larger N*/Dor T varies almost 
linearly with N*. We will see shortly that this variation 
is characteristic of a GausSian absorption and emission 
spectrum. For very large NSP only phonons in the far 
wings of the spectrum can escape and at these frequencies 
the spectrum is almost Lorentzian. Thus as before we find 
in this regime that the decay time varies as T, <> Ne 2, 
Since experimentally we do not find the linear variation 
of T. with N* at intermediate densities we conclude that 
inhomogeneous broadening is unimportant for the 29 en 
phonon line. 

As a final simple calculation let us find the decay 
time T. for a pure Gaussian spectrum for large N*. In this 
case 


a(v)b = a(v)b eget (5.45) 
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(a) Dashed lines: Slow decay time To vs density 
N /D, for various amounts of inhomogeneous 

broadening. Solid line: pure Lorentzian, Av ,=0. 
(b) Pure Gaussian. (Note the different horizon- 
tal@sca le Del Sederi nedminm Elm 2-4 6 O)h-eMmeIr Or 


G 
all curves T,¢/b =1/200. 


fee) 


where 


a(v_) = (5.46a) 


and 


AVG = 
De = moe /Th 5 Divas: CoO) 


Sopbq send 5B) —(5. S83) abecome 


a Ina) KOYAnKal(ve )b) 
a) - when r C5 247) 
0 Atyiee Yin (a(v_)b) 
yi 2 
Bee | n 2 1 i} 
soll is SS) hes dne Ao ee ° 7 
- V7 Th MO ye Velo ou) 


O (5.48) 
and 


saul: ra(v)b -¥in(a(v_)b) : (5.49) 


EX CecOUmEO@ECOluS tants Om OLOdeGraunLtye ChiS 1S ust the result 
found by Holstein? and Veklenko/! for Gaussian spectra. As 
in the Lorentzian case, Eq. (5.48) clearly shows that T. 
depends on the area under the wings of the emission spec- 
trum where phonons escape ballistically. Again spatiai 
diffusion is too slow and is thus irrelevant. Note that 
the radical varies so slowly that T. varies almost linearly 


with N* as we asserted above in the discussion of the Voigt 


eed) 


- 


— 
ply mo Sineqsk 
= 


dey 


profile. The exact result for all densities for a pure 


Gaussian is shown in Fig. 5.6b. 


oS oeebDascussion 

We have solved the set of coupled rate equation (5.1) 
describing the time evolution of a system of two-level 
atoms and phonons resonant with them, paying particular 
attention to the conditions obtaining during the final slow 
decay, i.e. during the phonon bottleneck. 

One of our most important algebraic results is 
Roc Cone 4 maw LeimoGLVvera simple algorithms for 
the slow decay time T aSwambhunctd Ol wOlmetne CONnCentration 
N* of scatterers for any lineshape. Notice that for large 
N* these equations reduce to the formal solution (1.8) of 
Holstein's integral equation. Because for large N* the 
phonons from the center of the lineshape are effectively 
Spatially trapped, T depends on the area under the wings 
of the lineshape where the phonons can escape ballistically. 
The system decays via ‘spectral diffusion' with phonons 
from the center of the line being absorbed, their energy 
redistributed over the entire lineshape on emission, and 
finally phonons from the wings (obeying a(v)-b< 1) escap- 
ing ballistically. ~Hence the system decay rate 1/T. equals 
the atomic decay rate 1/T, Cimesmthemrraction oLe the 
emitted phonons that escape ballistically. 

During the slow decay of the system a quasi- 


equilibrium obtains-", characterized by the two-level 


; 5 ee 
ef ice pyei! 


atoms and all the phonon modes obeying a(v)b>1 reaching 


SU cOoLmmPpOpULatLonsOLeavalue Na (tc) /N*. = thus lt, fOr 


2 
example, a(v) is Lorentzian the phonon population inside 
the excitation region is roughly rectangular with width 
Av = Av ova(vo)b <7 Ne? BYE Sle Gia Bisiig delat he 

Let us now specialize our discussion to the 2A¢+E 

transition of the cr?* ion in ruby. As evidenced by the 
experimentally observed dependence ee N* 2 LOnMlLarge N~, 
the wings of the 2A +E line are Lorentzian. There is 
insufficient inhomogeneous broadening to a Voigt profile 
to result in a linear dependence of T. Onan *<. 

Our model applied to the phonon bottleneck in ruby 

(with Lorentzian intrinsic emission spectrum) reproduces 
the following experimental features: 

(1) the rapid initial decay of the system in time 
tar The 

(2) the establishment of a bottleneck with a quasi- 
equilibrium excited atom population N, (eq) /N* 
Ne 2, 

(3) the establishment of a broadened (Av 2N**Av 0), 
flattened phonon population which has been 
observed by Dijkhuis et apet? in CW Zeeman 
splitting experiments, and 


(4) the subsequent decay of the quasi-equilibrium in 


1 
times He ee Nwageton tlatge Na. 


oe 
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om 
Wqyag ore? iau mi 
aa oie 
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To obtain quantitative agreement with the experi- 


2 = 
Ngee WeemustecnhOose tie 2A state 


ments byakbauli eta lz. 
lifetime Ty * 20 ns, if we take their estimates of N* and 
b at face value and take the density of phonon states 
D(v) accoLcdingge_O) EG wal(o.4) 7. LE we assume that -N* is 
actually underestimated by a factor of 3 (which is within 
the estimated uncertainty“+) and that the Debye approxima- 
ELOnMO D(v) is overestimated by a factor of 4, thena 
lifetime T,=9 ns gives good agreement with experiment. 


In any case our T, is larger than the values in the litera- 


1 
26 ; ike) s 
ture: Lengfellner et al. and Geschwind et al. , T,=-4 
ns; Rives and Neen T, =1.1 ety Wekals, Gbbaebls tah Age, 
T, = 2.8 ns; which we note already differ by a factor of 


seven and which with the exception of Rives and Meltzer, 
are all only indirect estimates. (Sees Ret. =13> for a 
critique of these values.) 

Rives and Meltzer?? treat the phonon bottleneck also 
via a rate equation approach. This model however pre- 
cludes spectral or spatial diffusion and assumes that the 
system decays via phonon decay. They evaluated only the 
early time behavior of the system and used the concentra- 
tion N” and the anharmonic phonon decay time on as free 
parameters. Their best fit was obtained with age, ae ns 
which is far shorter than estimates of 2 us previously 
sae ed Ne Furthermore their model cannot repro- 


duce the experimentally found functional dependence of the 
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quasi-equilibrium population N, (eq) /N* or the decay time 
T. aSmecuUncCELONS OG N*. 

A second possible resolution of the lifetime discre- 
pancy is the following: In a phenomenological theory such 


as ours, T., is the effective time in which, due to the 


au 
decay of an atom, a free phonon is generated. Now for 
Ne sD 2 we? aa", the concentration of scatterers is so 


large that a second atom has already reflected back a 
phonon before the first atom has finished emitting that 
phonon. Thus a cooperative effect can exist between the 


decaying atom and its neighbors, resulting in an effective 


decay time nee 2 Ty: This effect will be explored in 


Ghapter 6. 


CHAPTER 6 


PIEGD THEORETIC DESCRIPTION OF THE PHONON BOTTLENECK (IT) 


In Chapter 5 we solved a set of coupled rate equa- 
tions for the time evolution of the electron-phonon system. 


These equations were derived from the quantum equations of 


nNotronm4 oa mandm (47 jmeunder the assumption that N, < 1 Oiae 


cm > so that He OPN, T, and T., being the time of inter- 


action and the time between interactions respectively of 


the resonant phonons with cr?* ions. We have already 


indicated that under the opposite conditions when T, (N") ae 


: ; ama te? : 
cooperative effects between a decaying cr? TODAnGE LES 


i 


neighboring unexcited ions may exist, resulting in an 
increased effective lifetime T,- To explore this and 
other guantum effects not included in the rate equations 


we wish to study in this chapter the full quantum equations 


CAR >a) meanicles (4. /.).. 


6.1 General Description 
Let us begin by writing down the Fourier transforms 
of the Heisenberg equations of motion (4.5a) and (4.7) for 


the annihilation operator by (t) of the phonon mode of wave- 


~ 


number k and the lowering operator S(t) of atom n, namely 


g ; -ik-r, 
BU) aoe Cn Vian (x) (6.1a) 
k x OK aa) n 
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In these equations we have dropped terms containing the 


phonon operator b, (0), under the assumption that initially 


~ 


at time t =0 only a single atom is excited and no phonons 
are present. We assume that there are N two-level atoms 
at positions Cylons er lye Recall that the first term on 
the right hand side of (6.1lb) describes the spontaneous 


decay of atom n and that the second term describes the 


influence on atom n of the excitations coming from other 


atoms n'. (To derive the rate equations we dropped this 
term.) . We can solve (6.1b) formally by considering it to 
be a matrix equation with indices n,n'. Define the NxN 


path matrix (all matrices in this thesis denoted by 


carats) 
it BGAN OC) ge BE eA) ya EO Ganley arly” 
P nt (%) = d Gre 2)) 
ny 0 WTior ian ee, 
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is (the Fourier transform of ) the amplitude that an ini- 


tially excited atom is still excited, and 
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is the amplitude that a phonon is emitted by atom n' and 


absorbed by atom n. We can write (6.1b) as 


N 
$7 (x) =i a2tOM (x) 57 (0) ee } Lae 
n'#1 


(x) 50, (scien (6.5) 


This matrix equation can be formally inverted to yield 


~ “ oes "atom = 
m, (sc) Se leap) qe (SO);S mee (OPN. (6.6) 
n> Noy ny 
Here 7 is the unit matrix 1.. = 6... 
1) 1) 


Eq.) (6.6) has the following interpretation: ice (0) is the 


a 


probability amplitude that atom n, is excited at time t=0. 


i, 
(1-B)-* “4 is the probability amplitude that an excitation 
oe ae 
leaves atom n, and reaches atom n.. poem is then the 


1 2 


amplitude that atom n. does not decay. The product of 


2 
these amplitudes gives (after a Fourier transform) the 
amplitude that atom no iemexchredmatatimesr t... 16d ..1(6,.6) 
can be given an interesting pictorial representation. 
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where we have made the binomial expansion 
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In other words the amplitude that atom Nn, becomes excited 


is given by the sum of the amplitudes over all possible 
paths by which the excitations can go from atom ny to No - 
ei (SGee) cinch (dO) Clavel (elnyeye Sy is the amplitude 
for an excitation wtomnaves exactly sipineeracti1Onsebertore 
sabeskiaigte) ks, Ghelehi Nos the emission from atom ny Counting sas 
thewtvrstlinceract iron. 


Notice that in (6.10) an arrow representing a phonon 


cannot begin and end on the same atom (without going 
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through at least one other atom first; cf. (6.2)). The 
reason is that those diagrams would represent self-inter- 
actions causing the process of spontaneous decay, and they 
have already been included in the amplitude APs CLOmausL 
orders. This can be seen by denoting the amplitude of a 


free uncoupled atom by 


atom ik 
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Tie) roligtel Sa leale i.e. to inverse-Fourier transform 
Cee Ol 00) fee MNO FO mula t 21 ales can be expressed as some 
complicated function of A#tOM(x) and S22 63) Thus 4, (x) 
must have a branch cut along the interval O<x<)- 
Assuming that the resonance frequency A lies well inside 
this interval we can ignore any poles along the real axis 
outside this interval and evaluate Sa (Emus ingethesce— 
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(6.15) is our fundamental equation giving the probability 


amplitude that atom n, is excited at time t given that 
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atom ny alone was excited at time t=0 in the electron- 


phonon system. Our program now is to evaluate the ampli- 


tude MI and the amplitude matrix nee WiecheOoCcCurmeL netic 


integrand of (65-15). 
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From Chapter 3 we know that the amplitude A (xis) 


(where OSX 5 O57 6>+0") that the atom is excited is given by 
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We recall that the density of phonon modes D(w)dw per unit 
volume in frequency interval dw is given in one dimension 
by 


Diw)dw = aw (6.17c) 
@ 02 
D 
(2 is the length of the crystal), and in 3 dimensions by 


D(w) dw = } Di) (4) dw (6.174) 
=) 


where the density of branch s is (see Eq. (5.4)) 
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The Phonon Amplitude Matrix APA (3) 
Let us now evaluate the amplitude A, nt (xtie) 
Ud 
(O<x< a); 60") that a phonon emitted by atom n' is 


absorbed by atom n. We first do the 1l-dimensional case. 


1-Dimension: 
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where Si(x) and Ci(x) are the sine and cosine integrals>> 
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In the limits of small and large separations irae | we 
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3-Dimensions: 
In the 3-dimensional case we use the geometry shown 


below where ai a is along the z-axis: 
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In the continuum limit we replace 
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is the density per com> of phonon modes of type s in the 
interval dw, and a) is the speed of sound of mode type s. 


The phonon amplitude (6.4) becomes, after integrating over 


9 and 9, 
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An nt (xtie) = L ain ayeeene) F (6.25) 
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where the amplitude for phonons of type is is 
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limit of small and large separations respectively this 


amplitude becomes 
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Let us make the usual approximation that only fre- 
quencies near resonance are important so that we may 
replace x>A everywhere except in the phases in the phonon 
amplitudes Ae Conve nme Va Onc eeanOme oO. Zee tObesL—— andu3— 
dimensions respectively. Then we get for small separations 
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in both one and three dimensions, and for large separations 
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iirc ree dimensions. me LieEos. 9(On20)—(6.30)), A” and [ are 


given by (6.17a-e). 
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Notice the following features of (6.28)-(6.30). We first 


suppose that [eee ry >> c/A =A (A) /(27), 1.e. that atoms n 
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and ' are separated by a distance large compared to a 


resonance wavelength A(A). Then the phonon amplitude 
aio 

aph : sole ~r i 

An ni (tie) has the form of a plane wave we n on 
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in i-dimension and of a spherical wave 
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io /|\r iro in 3-dimensions, propagating from 


‘ CLS 
Dye 


i.e. that atoms n and n' are separated by a distance less 
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atom n' ton. Next suppose that Peta eo <aC/ Ai 


than = times the resonance wavelength. Then the ampli- 


tude qph 


a ni (%) Has the Giorm of a) virtual’ phonon'cloud 
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rather than a plane wave. This is evidenced by the fact 
“pn ; . Lf 
that AD ni (xtie) sesjustethe Shi ieenA ie Lorie of the 
J 
resonance frequency of a single atom undergoing spontan- 
eous decay. Thus at these short distances the two atoms 


nandn' sit in the same phonon cloud and we expect that 


the evolution of the decaying atom is greatly affected by 


the second atom. Indeed we will show presently that under 


these circumstances we regain Dicke's results on radiation 


trapping (see Section 2.2). 
i : atom : 
Having evaluated the amplitude A and the ampli- 
tude matrix apo in both one and three dimensions, we 
proceed and evaluate the matrix quantity (c= bisa come 


occurring in the integrand of (6.15). In one dimension 


we get 
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where the separation Pee =|xj-r, | of atom i and atom j is 
assumed large enough that the phonon amplitude has the 
form of a plane wave. Note the following features. For 
frequencies x far from the resonance frequency A' the off- 
diagonal elements of (ieeaP) are small and the binomial 
series expansion of (es jae will converge fairly rapidly. 
Physically this means that phonons far from resonance may 
escape from the system after interacting with only a few 
atoms. Various sneer have studied the simplest model 
of a single atom perturbing the lineshape emitted by an 


excited atom (the single particle approximation). This 


model gives a good description of the wings of the line- 


shape but breaks down near resonance. Exactly on resonance 
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the ij element of oes given by 
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with modulus one. Thus on resonance the binomial expansion 
of Gigs will not converge because all numbers of inter- 
actions become important? ’. The only recourse then is to 
TOOKmLO DEC cd GiaeSO Mil OnSmOl (O-3)))) ee ee SOLUCLONS. to all 
orders of perturbation theory, which we do in the next 


SecelonLtecnould peunoted that in three dimensions, the 


ij element (i#4j) of (1 -B(x+tie)) on resonance is given by 
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Although this matrix element is inversely proportional to 
the separation = of the atoms, and therefore is small for 
low densities of crea ions, it is compensated by the fact 
that the number of atoms ina shell of radius Poe about 


aon! Ss proportional to ae Thus also in three dimen- 


Sions we expect the expansion of ie 2 ee to diverge on 
resonance. 

Besides the divergence problem, the expansion of 
ak Sires does not conserve probability (i.e. the number 
of excitations or energy of the system) if the series is 


cut off at any finite order. For example cutting the 


series off at first order, we get 
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Pc icm Gu aaieeoru Cie Dime endil(6.16) for A 
(6.15) we find that the probability that atom j is excited 


at time t given that atom i was excited at time t=0 is 


given by 
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regardless of the number N-2 of other atoms in the system. 
But this is nonsensical because if there are for example 

7 atoms in the system under conditions identical to those 
of atom j then the probability that any particular one of 


them becomes excited must be proportional to WS 
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Evidently the higher order terms of the expansion of 
ii mye are essential in that they contribute probability 
amplitudes which interfere with the lowest order amplitude 


in such a way that overall probability is conserved. 


6.2 Exact Results 

We saw in Section 6.1 that a series expansion of 
oie in powers of P cut off at any finite order diver- 
ges as well as gives incorrect results. Thus in this 


section we attempt to find Gh yee Sxactly am elhisscansbe 


done for several special cases. 


Two-Atom System in One Dimension 
Consider the case of two two-level atoms, labelled 1 


Bude separated by a distance x. . = Sept | ee Tey Ay ime. 
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large compared to the resonance wavelength. In this case 


the phonon amplitude is a plane wave and we find that 
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we may also find it by expanding in powers of P to give 
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given that atom 1 was excited at time t=0, is given by 
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We can compare the results of the exact inversion and the 
series expansion inversion of (1-P) by looking at Figs. 
Opimpands6e2.5) fic, 0. lanandsb) respective hygshnow the quan— 


tities F, (x) and F, (x) (which are proportional to the 


Fourier transforms of the amplitudes S, (t) and S, (t) of 


the two atoms) with the matrix (1-py7t calculated to 


various orders. For example in Figs. 6.la, 1.3 the series 
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(a1 oes) TneSFourier transform of the probability 
amplitude that the initially excited atom 1 is 
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and second atoms are excited at time t, calculated from 


Picmcomrecnoondingequantits es, &. (x)rand F(x) Orel OS..0. la 


ll 
shekel Jel 


Notice the following features. Eq. (6.39) shows that 
Only even (odd) powers of P in the expansion of (l= Be: 
contribute to the probability amplitude for atom one (two). 
Fig. 6.2a shows that to zero™ order in the expansion, 
atom 1 decays exponentially and is never re-excited. To 
second order, atom 1 has a probability of being re-excited 


2ry 


once more at time t= aa , which is just the time it takes 


AeOHOnOHMtOmLbave le tOoratom 2 and back. “Keeping all’ orders 


in the expansion of Gomes yields the result that atom l 


: 2X12 
is repeatedly re-excited at times t= 


eal) eatin ely, 28) cen ers 
as the phonon is repeatedly reflected between the two 
atoms. After each transit the probability of the atom 
becoming excited is less because the phonon has a finite 
probability of being transmitted as well as reflected back. 


Notice that the single pulse in Fig. 6.2a, is higher 


2 
than the same pulse of Fig. 6.2a,- This is an example 
showing that the probability (of the first re-excitation) 


v1 1s cut off after second order. 


is incorrect when (ipeap) 
Evidently the higher order terms in the series reduce the 
probability sass wemsaldea tteor mono n. 

All features of Fig. 6.2 can be understood in terms 
GL the bourien trans: orms ot athe probabil lity, amplitudes 


shown in Fig. 6.1. Let us discuss the Fourier transforms 
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the transform is just a Lorentzian leading to simple expo- 
Monte laimdecay me lOmsecond order Hast (6.c07) and (6.39) 


show that the Lorentzian has an oscillation 


e (6.42) 


added to it. The envelope of the oscillation is largest 
On resonance and the maxima occur at frequencies satisfy- 
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Le. che ee peak in the spectrum occurs at a wavelength 
such that exactly n wavelengths fit in the distance 2r1>- 
In other words a peak in the spectrum occurs when a stand- 
ing wave is set up between atom 1 and 2, the height of the 
peak being the greatest for eee near resonance. 
Binally sncludingealiechesb1 ghemsnumbocr smo pe pnOonOnm re .lecs 
tions gives the spiked spectrum of Fig. 6.la,. Simple 
Fourier analysis gives the result that the spacing between 
the peaks in frequency space is inversely proportional to 


the period of the system. Specifically the re-excitation 
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The narrowness of the spikes in Fig. 6.la, yields the 
Tac oCenNunper Or Le-excitatron OL el tq. 6.2a,. 

An important case occurs when T, > 2rj,5/c, i.e. when 
the lifetime T) of the decaying atom is longer than the 


time it takes the emitted phonon to be reflected back to 
the atom. Then the Lorentzian is narrower than the spac- 
ing between the spikes, and the decay of the atom depends 
crucially on whether or not the Lorentzian overlaps with 
one of the spikes. In the case that it does not, the 
resulting lineshape is approximately the original Lorent- 
zian and the excited atom decays in a time t + Ty more 
or less unaffected by the second atom. However in the 
case that the Lorentzian does overlap with one of the 
Spikes, the resulting lineshape is much narrower than the 
Original Lorentzian. Physically this’ means that the reson- 
ant phonon has formed a strong standing wave between the 
two atoms and the two-atom system decays very Slowly in a 
time Tea 2 ila 

Let us consider whether or not this phenomenon can 
OCCULT In OPtLCa! Lye DUMDCOMLUD me CEN Ss Ce DeCenOcCCds Liat 
since ruby is a three-dimensional system, the amplitude 


for a phonon to be emitted by atom 1 and absorbed by atom 2 


is given by 
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This is the same as the one-dimensional amplitude 
(C2) except schatminmtnis case the: amplitude 


decreases with the atom separation r Therefore for a 


2s 
resonance to occur between two atoms in three dimensions 

we have in addition to the former standing wave conditions 
that Kyo =n Vi2maned Ti 


SEC, Ah eC sethatuthe separation r 


<T,c/2, the new condition that 


ae is not much larger 


ied iF 2 
than a resonance wavelength. The wavelength of phonons 


resonant between the 2A and E levels in ruby is ro ~ 10 “om. 


Assuming that the average distance between neighboring crt 


POnsmatbwasGensttyaNpeis ec tau. > = yer? 


Aen , we see that ata 


=}_, 1.e. nearest neighbors 


pemeouLyaN es l0ascman <r > 
Tel Ws fe) 


are on the average a resonance wavelength apart. Thus at 
this density we expect many Cran ion pairs to satisfy the 
standing wave conditions. It should be pointed out that a 
proper calculation would involve studying the radial dis- 
Eribution function torsexcited Cra ergs? givay veibleyy7-, — “UMley=s 


Gadialedastributionmeunct1 on no(r ) is proportional to 


eZ 


the probability that given a Cra ion is located at some 


1’ a second ion is located a distance Yi from 


that ion. This function would depend on how the crystal 


DOsi tone 


was grown, etc. In general no (rj5) is a maximum near 
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Loe on a >. If we assume, as is the case in liquids? ’, 


that this maximum is rather broad, then even at densities 


lower than N* = ihe com> we expect the standing wave con- 


dition to be satisfied for many ion pairs. Although more 


study is required, we conclude that this phenomenon could 


account for the long effective lifetime Tei noted in the 


rate equation analysis of Chapter 5. 


N Atoms, All at the Same Position 


Another case that can be solved exactly is the system 


of N distinguishable atoms all at the same position. In 
practice this would occur if the atoms were more than a 
deBroglie wavelength apart (i.e. distinguishable) but less 
than a resonance radiation wavelength apart. This is the 
system studied by Dicke?? (see Section 2.2). It is easy 


to verify that for this case the inverse of the matrix 


(lhe Gee aoa (ome) 
jie ety aS 
where 
a=x- ated ; (6. 46b) 
and 
pb=a-at+it ; (6.460) 


is the matrix 
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Ube 4a) 
where 
N-1 ale 
N 
c= i+ 4 = ; (6.47b) 
<i en A ih te ) 
and 
_i 5 
N 
dS sae 2 a (6.47) 
cept CASS i te) 


Plugging) (6.4 /a,b,c) into (6.15) we find that the probabi- 
lity that atom n is excited at time t given that atom 1 


was excited at time t=0 is 


it Lyre + - 8 
Niet) 605 Sn ft) 53 (4) | (5° 
: : 2 
= eee ee a Nt NE C/ 2 
(Sy 4 nN? ER e | 
( =a 
a oe P FOr tes osUar, fLOTsM a= elk; 
ia, 
= Ce ; idoke ae Su sehen chean I, (6.48) 
= : for t > »~for all N-l other’ atoms. 
N 


(6.48) shows that the probability of non-decay of atom l 
decreases initially at the rate T of a single isolated 


atom but soon becomes fixed at the non-zero value 


a (eco ar 
Nj (t ) ( 


probability 1/N* of becoming excited. Hence the probabi- 


———) . Similarly the other N-l atoms have a 


lity that a phonon is eventually emitted is only Say ee 
This is in complete agreement with Dicke's results on ra- 
diation trapping which we presented in Section 2.2 

Let us look at the phonon spectrum that is emitted. 
Plucoundacon djmandm ono) 2ncCO Ne Fourier transtorm for— 


mula (3.3) we can write 


Saec0)) . ea 3 A la 
b, (t) = 3-5 mae THe © pees aak ees ee F 
J ie ageuk n,1 
£ 
(6.49) 


where © is the contour (3.3a) around the lower half- 
plane. In this equation, since all atoms are at the same 


DOsLeLonyewe find that 


. res atom ORT ce ah JE 

e A (xtic) (1l—P (xtic)) So Gime re oe eee 
a= ee ee) 

(G75.0) 

Thus doing the contour integration we find that 

S_(0) g -iw,t -i(A-N(A-A')t -NIt/2 
ih k 
b(t) = iT (e -e e jt 
ee A+ N(A- At +—-) 

(6,025 45) 


The total number of phonons in the system at time t is 


given by 
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2 1 + 1 
Noy (t) = Ls (5), 1b), Ce) b, (t) GLE 
0) 
: ( 
é By Weis aa te sae _ 5 NIM? 
" SSN eer.) 
x cos (w= A +N(A-4")) th (6 .52a) 


By either expanding the quantities in braces for t=0 or 


by integrating over-@ we can write (6-.52a) in the following 


forms 
(7) aS 
| dw =— ie eOne fe St 0) C625 25) 
f T 
Napa Set l= 
l -NTt ts “Pemcor et a~t0 G6r. 2c) 
Wy eee Pate 
N r LO Gants i £3 


Eom (G7 2b) )meshowsmthateror.stimest = 0 vasunbformuspectrumiof 
phonons is emitted in agreement with the uncertainty rela- 
Elon 0) 1tee-al eo sob) smradentical to the result for 

the emission spectrum of a single atom for time t=0 (see 
GrEsScussi OneOLeLd an 41L3))).0ee Oa 2C mSOOWSE tate lOm erin? be 
but small times the rate of phonon emission is I, again the 
same result as for a single isolated atom. Finally (6.52a 
and d) together show that after a long time the N atom 
system has emitted a phonon with probability 1/N, and that 
the phonon emission spectrum is Lorentzian with width and 


shift N times that of the Lorentzian emission spectrum of 
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a Single isolated atom. This is an extreme example of the 
well-known phenomenon of pressure broadening, i.e. the 
broadening of the emission spectrum of an atom due to the 


perturbations of nearby atoms. 


N Near Atoms 

Having obtained the exact solution to the problem of 
N atoms at the same position, we can do perturbation theory 
around that result and consider the case where the N atoms 
are only near each other. (By near we mean that BERS Cy 
c/A where MERI yo! is the maximum separation of any two 
atoms.) 

From a mathematical point of view the radiation 
trapping that occurred in the former exact case is due to 
the pole at x=A in the probability amplitudes (6.47b and 
c). When the atoms are separated slightly we expect this 


pole to move slightly off the real axis to a new position 


slow 


2 
Slowly leaks away at the rate Im(T 


x= - with the consequence that the trapped radiation 


i. 


slow 


Let us begin by defining the matrices 


n~ 


Ih nw as 


1, (*) - atone = P(x) ; (6m530 
A (x) 
and 
vas ic Ve x 
y(x) = —=— (1 - P(x)) . (6.54) 
OPEL ices 


ie and ) are respectively the zero-order matrix (i.e. the 
matrix of the former problem defined in (6.46)) and the new 


perturbed matrix. The difference between these matrices 


s 4 
r ; - at YR ~ = _= ¢ £* 
‘ A a ** - @ awe & 7 


‘ * a 4 - ae a 
4 39OC "O77 28 1 $i32 Oo PAANSOSOT 
« ; ¥ —_— = 


~ = re Ta 
(L655 Pf 


defines the matrix 


e(x) = ) (x) pales) peecanre (Bo (x) = P(x) ) 
x 
= APP (x) - APP (x) (ess) 


, can be inverted by using the fact that 


a eo ee (6256) 


It is convenient to define the matrix E, such that E,,-1 
ih 


io Cease, | anast ne quantity Bos et AE ec Then we may 
write 
Van (any SE erin ale 2 (xtie ) yee 
Saul. ape B 3 
~x-A ie x -A +NB6 B) ra 
and Eq. (6.56) becomes 
ot tic) = ar OMixtie) (1 - B(xtie))” 
ek 
eee Aas: aoe te pape 
Tene GG PTI x= A) feet, 
where 
v= = Qe - a cE Ge OC (6.59) 


Ss 
Thesshite Of thes pOleut Lone x.=(ieLO x=A-+- is described 


by the formula 


eZ 


z des bese’ Z ro 72 s,2 
emer Ss Caen At (1 Spank + mOEG ae aye)! he (6.60) 
2 


Comparing (6.58) to (6.60) we find after some algebra that 


[erie AG eye Pp ttie)) ae 
ag) ij 
= 
ij N i: alee NE (ena 
Ss eT ( ; 
ves Re Meh oA eS) 
xXx— At an 


where the shift of the pole in the ij matrix element is 


given by 
il di 
po, HE N dates) +5 ) fxg 
ei) Bese aticke yor SKN ake (6262) 
2 eee ote 
ij oN 


with a) aa evaluated at x= A. 


Two Near Atoms 

Consider the simplest non-trivial example, namely 
that of two near atoms in one dimension, with atom 1 ini- 
Evel yee <e) todpanag atolmle2 sinitiallyounexcited.,=iIn this 


case it is easy to show that 


ie = 

ae os Re (€, 5) = if (4S) r (6.63) 

rj54 
(for 7c << 1, where Yi is the separation of the atoms) 
and that 

ie. by 7 

S, ets eet 2 

Tae = 2E45 214--( 7G ) P (6.64) 
for all ij. (We ignore the real part of the shift Tia) 


Conmbin nga omo4mwl the ovo banda Ondo jeewer find that the 
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Prana owuloveithdtratOMen,— il — i, 2, luseexcrted at time t, 


given that atom 1 was excited at time t=0, is given by 


ay 
D : =a, —J 2 
4 1 ca ieee | 
N = atte ~ - i 
nae) a | dxe ae: ic) Jai (x+ie) || 
0 
Meee t/2 2 
={(g. -2) gidt, slow ph gi(a-2(a-a))t cre| 
ee 2 
(6 .65a) 
The slow decay rate for both atoms is 
ion BAVA. 9 
a 12 
ieee = 2T ( Fc ) : (6.65b) 


BNeses probabilities sare shown, in Fig. 6.3a and b for two 
different separations Ti of the atoms. Note that the 


decay rate See, increases as the separation increases. 
The inset of Fig. 6.3(a) shows the Fourier transform of 


the probability amplitude for atom n, namely the quantity 


Vay (xrie ) = [oo Getic) 


ITT (67.610) 


We see that the former pole at x=A has become a narrow 


resonance. 


N Randomly Distributed Near Atoms 

Suppose there are a large number N of two-level 
atoms randomly distributed in a one-dimensional box of 
Vengthel, awitheonlyeacolmel mind tidis Vaexeci1 ted, a lulls cacy, 
to show that the average square distance ane between any 


two atoms is given by 
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(a) 
ELG=so woe a) nes probab iiiia es Ny (t) and No (t) that atoms l 
and 2 are excited at time t. System parameters: A = 
“401 T= .20p, and atom separation X72 =C/Wp: 
(b) Same as (a) with atom separation Ly 9 = C/U: 
The inset of (a) shows the Fourier transform of the 


probability amplitude for the two atoms. 
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<Xj5° = [A/6ee (Oreo) 


thus frome (6 763) 


2 
See Gk 
<Im(e),)> = G (50) : (6n60) 
If we now assume that 
<o > ; FOLs Leas) 
E4505 12 (eneg) 
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Formed eva. er NuSsmENeDLObabLIt ty ethatyatom 1 (is excited 


atetime tis 


ly ie 
Ny (t) Slow (ome 
where 
ee LAs 
ow 1 3 (So) : (6 .71b) 


The probability that any of the other N-1l atoms become 
excited is negligible for large N. Let us now find the 
phonon spectrum that is emitted in this case. The density 


(w, )dw, of phonons emitted into frequency interval da, 


eh! Keeeek 
is given by 


‘ 1 + il 
Dp (yr t) doy, = D(a, )2< (9) 4 [By (E) By CE) | (9) > : (6.72) 


where by (t) LSeOlveleby 6 Gn 49))) pelo uta nd b, (t) we must 


evaluate the quantity (see Eq. (6.50)) 
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In the previous case where all atoms were at the same 
position, the contribution from the first term in braces 
vanished identically. However in the present case with N 
large, this term gives the dominant contribution. Thus 
neglecting the second term and using the expression (6.70) 
for Tis we find that (6.73) becomes 

San ; (6.74) 


x - Atl _Slow 


where C is an undetermined constant. Plugging (6.74) into 


(6249)ewe find that 


i ~lo,t 
CgS, (0) e 
[Ojy (RESON ° (ores 5:) 
k ey 
ep, i Se 
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The condition that a phonon is emitted with certainty as 


t>° determines the constant Ic}? EOEDe 


r 
Kink SE , (omen 


The phonon spectrum emitted after long times is thus 


ls low 


Xai; 
aL genie rete = : (Gi) 
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Let us paraphrase these results. If a single atom of a 
collection of N atoms randomly distributed along a length 
L < c/A is initially excited, then that atom decays at the 


rate 


"siiiie aa SURE rae 


emitting a Lorentzian phonon spectrum (6.77) of width 


r This decay is indistinguishable from the decay of 


etiae 
, & 


slow’ 


a Single isolated atom whose effective lifetime is T 


Ak Meet smincecoouingueto more that thes rate 


is 
slow low 


independent of the number of atoms N (for large N) in the 
system. If we imagine that the length L of the system is 
Musee Nemrange Ore tne svittual, phonon) Cloudssurrounding the 
excited atom (i.e. the range for which the near-field 


form of the phonon amplitude (6.28) is valid); namely 
ky SRY GGri2)) 
then we find that 


mae 
etivens ea |) : Hane) 


This result holds whenever the number N of atoms in the 
range (6.79) is large compared to unity. 


Assuming that these results hold @also im the) three- 


6 & 
dimensional ruby system where L = 13 _amesec ANG H cm, we 
2a Oma Zz 34. 
expect (6.80) to obtain for densities of excited Cr ions 


n* > Wo x Toe cone Unfortunately these densities are 


beyond the range of the present phonon bottleneck 
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experiments in ruby. To our knowledge this phenomenon 
has not been looked for or observed in optical radiation 


trapping experiments. 
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CHAPTER 7 
CONCLUSION 


We have studied the phenomenon of resonance radiation 
trapping by a system of two-level states paying particular 
attention to the example of the bottlenecking of 29 cm + 
phonons by excited Coat LOnS sungoptuca bivy tpumpederuby 7, In 
Chapter 2 we found the field theoretic form of the electron- 
phonon interaction in ruby, and in Chapter 3 we studied the 
time evolution of a single two-level state coupled to a 
phonon field. With this experience we proceeded to solve 
the problem of the time evolution of a system of many atoms 
eoupled to) a phonon Greld by deriving in Chapter s4)-a).set:of 
rate equations describing the system. These equations were 
shown to be valid for small numbers of excitations and at 
low density of two-level states. The former restriction 
allowed us to neglect stimulated phonon emission and the 
latter to neglect the quantum interference effects due to 
phonon reflections by nearby atoms. In other words an 
electron-phonon interaction could be viewed as a two-step 
process with the phonon emission by an atom independent of 
the preceding phonon absorption. With these limitations in 
mind we solved the rate equations in Chapter 5 and compared 
the results with recent experiments on the phonon bottleneck 
in ruby. The rate equations were able to reproduce the 


following experimental features: 
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(1) the rapid decay of the system in time t x T, from 
its non-equilibrium initial conditions. 

(2) the establishment of a bottleneck with a quasi- 
equilibrium excited atom population N, (eq) /N* es nen and 

(3) the subsequent decay of the quasi-equilibrium at 
the rate 1/Ttoo 1/T, xf where f is the fraction of the emit- 
ted phonon spectrum which is not resonantly reabsorbed. In 
the case of ruby where the emission spectrum is Lorentzian we 
had £ = 1/va(vo)b. 

We saw that the ne? behavior of the quasi-equilibrium 
excited atom population was a consequence of the fact that 
all phonon modes in the frequency band satisfying the condi- 
tion a(v)b>1 came into equilibrium with (i.e. reached the 
same population as) the excited atoms, and that this band 
had a width Av VYatv,)b 2 N*?, 

The appendix shows that this broadened, flattened 
spectrum existing inside the system is the same as the 
Spectrum that would observed emerging from the system. 

To obtain quantitative agreement with the experiments 


PA Ny A 


DVebaltligeetrals we had to choose the 2A state lifetime 


T) ~ 20 ns in the rate equations, which is longer by an 
order of magnitude, than most previous estimates of T,- 
In Chapter 6 we showed that this large lifetime could be due 

to a quantum interference effect in which a standing wave or 


resonance is formed between the decaying atom and a neigh- 


boring unexcited atom. The duration of this resonance then 
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the otherwise classical rate equations. 


enters’ as the effective Litetime T of an excited atom in 
In Chapter 6 we found another guantum effect which 
occurs at densities so high that there are several other 
atoms within a resonance radiation wavelength, i.e. within 
PRS) ees 


eheevirtual phononeclouds ofsansexcited atom (N~ = 10 cm 


in ruby). In this regime the radiation, which is perfectly 


SEE 


trapped in Dicke's model, slowly leaks away in a time T) 


2 Ty) independent of the atom density. 

Some of our assertions could be tested experimentally. 
For example the radiation leaking result of the preceding 
paragraph could be checked in impurity-doped crystals such 
as sapphire which have a longer resonance wavelength than 
ruby, and for which the effect should therefore appear ata 
lower density N*¥. Also no experiments to date have examined 
the dependence of T OnmehemscezeeblOmitierexci tation region. 
We predict a b? dependence where a theory due to Pauli et a1 29 


predicts that To yee and purely spatial diffusion implies 


that Beet: 


Os 
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APPENDIX 


THE EMITTED PHONON SPECTRUM 


Until now we have not discussed the spectrum of 
phonons that emerges from the excited region of the ruby 
crystal. This spectrum can be derived from a simple 
phenomenological argument. Suppose that we have a uniform 
distribution of two-level atoms in the upper energy level 
along a one-dimensional system of length b, as shown in 
Fig. A.l(a). Then each infinitesimal region of length dx 
about x, 0<x<b, of the system emits a spectrum A a (v)dx 
of phonons toward the observer located at x=0, where A is 
some constant and a(v) is the frequency dependent absorp- 
tion coefficient (which is proportional to the emission 
spectrum). In travelling to the observer this emission is 


attenuated by the factor ema 


Integrating over the 
entire system we find the spectrum P(yv) arriving at x=0 
from all parts of the system, namely 


b 
P(v) =Aa(v) | cuore ausulcur ani (weet) Dye ca) 
0 
Ef-o.(vjb eeelyirves ifethersystemeis So transparent, 
then P(v) =Aba(v), which is just the unattenuated emission 
Spectrum. sOnethevothewahands liao.) mo cL pc aL oe tie 
system is opaque then, as shown in Fig. A.l(c), P(v) =A= 


constant over a broad band of frequencies around resonance 
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observer | Te 4) 
x x+dx b 


(a) 


ae : AV, —> 


(tor, LOoreneZziangnc cy ).) 


(el) 


Fig.A.1l (a) Geometry of the one-dimensional radiation 
(b) The absorption 


trapping system of size b. 


coefficient a(v). (c) The phonon spectrum P(v) 
het eee 


emerging from the system in case a (vo 


PFI 


Satisfying a(v)b>1, and P(v) =Aba(v) in the far wings 
where a(v)b<1. In case a(v) is a Lorentzian of width Avo 
then the emerging phonon spectrum is broadened to width 
Valv yb Avy and flattened similar to the phonon spectrum 
inside the excitation region shown in Fig. 5.4, except for 
the absence of the dip near resonance. This self-reversal 
Gfgeuhetd Stributvonein rig. ©5 -4iis due to the*vanishing of 
the density of Cran ions in the 2A state at the walls of 
the container. We can reproduce this by allowing the cons- 


tant A in Eq. (A.1) to have a spatial dependence such as 


that of the system of Chapter 5, namely (see Eq. (5.12)) 


A(x) = J. (=-|x-s]|) . (A. 2) 


Thus we find that in the steady state the phonon spectrum 
emerging from the excitation region is the same as that 
inside the excitation region. 

Some features of the emitted phonon spectrum can be 
found via quantum field theoretic methods. Recall from 
Ghaple peom(EC sie 6 ala)sand 6.0 peetnatethe probabiliry 
amplitude that a phonon of wavenumber k exists in the sys- 
tem is given by 

N -ik-r 


B(x) ==2— Je ~ "gi (x) , (A.3) 


aa 1 n=1 
where the amplitude that atom n is excited is given by 


=i shee 


Bo (x) = (1 - P(x)) A 2S) Ae) - (A. 4) 
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ee 


We have assumed that in the initial state je eh es 

atom m is excited and that no phonons are present. From 
our work of Chapter 6 we know that whenever the two-level 
atoms are not all at the same position, then 4 | (x) has a 


branch, cut along thevreal axis in the interval 0<x<w4 


D 
and no poles. Hence its Fourier transform S(t) decays in 
a finite time. Since we are interested in the spectrum of 


phonons existing in the infinite time limit (when all the 
phonons have escaped from the excitation region), we may 


2 alge | b alt...) by simply evaluating B, (x) at the pole x=w 


k k’ 


Nt Sane anes ws atom - 
Dett+>) = lim ) ~e Clam ti shoes es (a tie)S (0) . 
e>0 n=l n,m a 


~ 


(A.5) 


The number of phonons eventually emerging from the excita- 


tion region is given by 


(to i>. (A.6) 


Bie i ae 
Non pe, Pit ibe 


Let us now specialize to a one dimensional continuum 
system. Then (A.6) can be written 
D 


W 

Beh 

hss = | Megat te! dw, A (AS 7) 
0 


where 


on au\ted 
% bs . >) = ms Sa a J > << . a 
| | +. %) tn mateye 


a 
t 


iw 
N k 
Lee Saag ais we ; ~1 atom 
Doh (4) 1 = ep ae ley, es (L-P(w,tie)) A (a tie) 
10 
| Nie ss fn Aven -l atom 
+} ) e (l-P(w,tie)) A (w, tie) A ARS) 
| k n k 
n=1 m 


Here the first/second term in braces represents phonons 


emitted to the left/right. Let us arrange the atoms in 


the order shown below 


(A.9) 


and assume that the distance Eye between atoms i and j is 
large, L.e. Eye >> c/A. Consider the spectrum of phonons 
emitted to the right by an interior atom m. Note that we 
may formally invert the matrix Gite Co Pile) meOCCUrEI NG ein 


k 
(A.8), yielding 


(A.10) 
(w, tie) ) 


where Cs (w,tie) is the cofactor matrix formed by multiply- 


ing eagye 2 times the matrix formed by deleting the qe row 
and the aan ere biiigy ese (1- P(w,+ie)). Thus we can write 


the second term of ({(A.S)) in the form 
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<1) 2o\pmuyloo ~~ 


iw 
atom N --=r_ 2 
Nop (@) do| = 5—|* ane) Cc MG (wtie) 
emitted to det (1-P (wtic)) n=1 
chem elon t 
by atom m (A.1la) 
atom , 2 2 
“Fa ee | PICSEMGNGtt eli eye yo(Aalib) 
det(1-P(Wwtic) ) 
where 
det C! (wtic) = th see 


NXN matrix 
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where 


det C" (wtie) = 


inn - aimee 9 XT. Mate Ky 


(A.13) 
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and 
a(w) = ES (A.14) 
' abel 
(oe \ ae 
2 
Eq. (A.lla) has the following interpretation: 
At time t=0 atom m is excited. 
Sim@liny Api , 
A is the probability amplitude for the non-decay 
Opethasmea tom. 
WAG. <8) ode (le pyre is the sum of the amplitudes 


over all possible paths (traversed all numbers of times) 
for a phonon to be emitted and to eventually return to the 
Crloiioll veces CLtEecracomeam. 

ee is then the amplitude for a phonon to be emitted 
by atom m and absorbed by atom n, by traversing a simple 
non-repeating path from m to n. The product of the above 
amplitudes gives the amplitude that atom n is excited. 
Finally the amplitudes for phonons to be emitted by any of 
the N atoms add in phase to give the amplitude for the 
existence of a phonon. 

Rose (hai D)jeand (Asl2)mashow thatpby (factorizing 
det C' a different interpretation is possible: 

detsC" is the amplitudes for) agphonon)toyrnteract 
with the m-l atoms to the left of atom m and to eventually 


return to atom m. 


The factor 
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F(a) = [t-a[*%™ = wis) | (A.15) 
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is then the factor by which these phonons are attenuated 
in travelling the distance from atom m to the right to 
acomeoNns. 

To date it has proved impossible to evaluate 
det(1l-P(wtic)) and det C"(wtie) in closed form for all w. 
However if we approximate these quantities by their values 


far from resonance, namely 


det(l - P(we~tice)) = 1, 
CAREEG)) 
det C"(wremtie) = 1 r 


then Eq. (A.1lb) becomes 


ST it 
n_, (w) dw at i Se Ls) 
ue m;right Ata (a —A')* + (5) ot 


The attenuating factor F fo) and the spectrum 


(w) emitted to the right by atom m are sketched 


eh 


in Figs. A.2(a,b) respectively. We see that the presence 


m; right 


of N-m atoms between atom m and the surface of the excita- 
tion region results in a 'hole' of width=IvV/N-m in the 
center of the spectrum emitted by atom m. 

If we assume that the one-dimensional system is 
uniformly excited, i.e. that each of the N atoms has a 
probability 1/N of being initially excited, then the 


averaged emission spectrum is given by 
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(a) The factor ERS) by which phonons of frequency w 


are attenuated in travelling from atom m past N-m 
atoms to the right-hand boundary of the system (see 
HGz\As15);) ab) B Thesphononmspect rum De eect 
emerging from the right-hand boundary of the system 
given that atom m was initially excited. (The light 
line is the unattenuated Lorentzian emitted by atom 
m) (see Eq. (A.17)). (c) The phonon spectrum ee 
emerging from the system assuming that each two-level 
atom has probability 1/N of being initially excited 


(seeuhoq. (A. 13))". 
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We see that the spectrum <M y (0) >s which is sketched in 
Hope Cu metsmer lateandenas widths N I. 

In this one-dimensional system a resonant phonon 
undergoes exactly N absorptions in traversing the system. 
Thus it approximates a three-dimensional system in which 
the absorption probability a(v )b=N. Thus this spectrum 
Gsesimilarn tomthe one of width Ya(vo)b T derived in (A.1) 


by the elementary argument. (However due to the rather 
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drastic: approximations (A.16) it is normalized incorrectly.) 


It is possible, by doing the calculation on a compu- 


ter, to find the spectrum eats 


excited one-dimensional system without using the approxi- 


mations for det(1-P) and det C" given in (A.16). The 


results for systems of N=1,3,5, and 8 two-level atoms are 


shown in Fig. A.3(a-d). Here a Monte-Carlo-type calcula- 


tion was used: Positions for the N atoms were chosen by 


a random number generator 128 times, each time the result- 


ing spectrum being found. The spectra of the 128 ensemble 


members were then averaged to produce the spectra of Fig. 


A.3. Also shown (light lines) are the approximate spectra 


(w)dw> emitted by a uniformly 
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(Ca) (b) 
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Fig. A.3(a-d) The exact phonon spectrum Sieh a emitted 
from a uniformly excited N atom system; N=1,3,5,8 
respectively. The light curves show Sag WRN accord- 
ing to the approximate formula (A.18). Note the dip 


in the spectra on resonance. 
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according to Eq. (A.18). We see that the approximation 
underestimates the population of phonons with frequencies 
within a distance YN 5 of resonance. The spikes in the 
Spectra are due to long-lived resonances between various 
atoms of some of the systems, and have not yet disappeared 
in the averaging process. It is interesting to note that 
the spectra have a hole at the resonance frequency. This 
Nees canmbemShownmeombe due tovthe fact) that in a one- 
dimensional system, phonons exactly on resonance are 
Bees cu packwardé from an unexcited two-level atom with 


100% certainty. Since this phenomenon occurs for all 


phonons except those emitted outwards by the two end atoms 


D 
1 and N, the hole has a height =< =x the hergntsor sthe 
Original Lorentzian. Since in three dimensions the reflec- 


tivity of resonant phonons is not 100%, this effect will 


not be observed there. 
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